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Abstract. We present global spectral energy distributions (SEDs) from 0.3 /im to 90 cm for a sample of seven 
low-metallicity blue compact dwarf galaxies (BCDs). In addition to data compiled from the literature, we report 
new SCUBA data for the galaxies in the sample, including the two most metal-poor star-forming galaxies known, 
I Zw 18 and SBS 0335—052. The standard starburst templates, M 82 and Arp 220, do not give a good approximation 
to any of the BCD SEDs in our sample. Several SEDs are instead characterized by: (i) approximately flat radio 
spectra suggesting dominant thermal processes; (ii) far- infrared (FIR) spectra which peak at or shortward of 
60 pm; (in) a signiflcant warm dust component and absence of Aromatic Features in Emission (or PAHs) in 
the mid-infrared (MIR). The SEDs of low-metallicity BCDs do not follow "canonical" scaling relations, and the 
"standard" star-formation indicators based on radio continuum, FIR, and MIR luminosities can be discrepant by 
factors of XL 10. We present new models of the dust SEDs from 5/im to 1cm, and derive infrared luminosities, 
dust distributions, temperatures, and masses. The observed SEDs and dust models are interpreted in terms of 
the active/passive regimes of star formation and trends with metallicity. Finally, we discuss the implications of 
our results on the z ^ 6 starburst populations which will be detected by forthcoming submm and radio facilities. 
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1. Introduction 

The physical processes which governed the formation of 
the first stars in the universe and the epoch during 
which it occurred are a main focus of modern cosmol- 
, ogy. Measuring the cosmic star formation rate (SFR) 
of the universe as a function of redshift or look-back 
time has been the subject of much effort, initia lly at 
ultraviolet wavelengths (e.g., iMadau et al.l Il996|) . and 
more recently in t he millimeter and ra dio regimes (e.g., 
iBlain et alJ Il999l: iHaarsma et~al] l2000|) . Unlike the ul- 
traviolet or optical spectral regions, emission at radio 
and (sub) millimeter wavelengths is unaffected by dust, 
which makes uncertain extinction corrections unneces- 
sary. To exploit this advantage, emission at virtually ev- 
ery wavelength from the mid-infrared (MIR) to the ra- 
dio has been used to trace star formation, both localh 
l)Kennicuttl ll998": Rousscl et al."200l'; 'Condon ct al."l99l| 
and at z 5, 4 ( Row-an-Robinson et al.i . 1997; .Blain et al. | 



- ICariin fc Yunl 119991: ICharv fc Elbad 120011 ICarretl 

2002tlYun fc Carillill200^ . Nevertheless, the conversion of 
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MIR, far-infrared (FIR), sub- mm or radio continuum flux 
to a global SFR depends on canonical scaling relations 
(radio/FIR, MIR/FIR), which may not be applicable to 
primordial star- forming systems at z > 5. 



Even locally, there are some indications that SFR is 
not adequately traced by FIR or radio lumin osity of low- 
luminosity or young sta r-forming galaxies jKlein et alJ 
Il99lt iRoussel et al.ll200^ . Indeed, in dwarf galaxies, the 
scaling relations appear to hold only because both IR an d 
radio emission underestimate the true SFR l|Bel]l l200^ . 
In these objects, particularly those of low metal abun- 
dance, the SFR estimated assuming the standard radio 
mix of thermal/non-thermal cont ributions can be under- 
estimatcd by a factor of 5 or mor e teobulnick y fc JohnsonI 

[T999.; .Beck et al. 2002; iHunt et al 2004a,). While rare in 

the nearby universe ( see iKunth fc OstlinI f2000|) . low- 
metallicity low-mass galaxies might be much more fre- 
quent at early times and high redshifts, give n the predic- 
tions of the hierarchica l merger models (e.g.. lBaugh et alJ 
Il998t ICole et al.l 12000*). Indeed, such objects may repre- 
sent the primordial "building blocks" - or "sub-galaxies" 
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l|Reeslll998|) - in hierarchical scenarios of galaxy forma- 
tion. 

To better understand the effects of metallicity and lu- 
minosity on the properties of star-forming regions and 
investigate possible spectral templates for high-z galax- 
ies, we have undertaken a study of the spectral en- 
ergy distributions (SEDs) of low-metallicity blue com- 
pact dwarfs (BCDs). Low-metallicity galaxies are in fact 
the most viable local laboratories with which to study 
primordial galaxy formation, because of their low mass 
and chemically unenriched interstellar medium (ISM). 
However, such objects cannot be compared to the mas- 
sive "SCUBA galaxies" and Lyman-break s ystems which 
are already chemically enriched at z ^, 3 jPettini et al.l 
l200lt iTecza et al.l l2004t IShaolev et alll2004() . At higher 
redshifts z ^ 6, truly primordial star formation must occur 
in very low or zero metallicity environments, and metal- 
poor BCDs are the only way to study such processes 
locally. Because IR-radio correlations are used to infer 
SFR for galaxies observed at high redshift, it is impor- 
tant to verify that all galaxies, even the most chemically 
unevolved, follow such correlations. Moreover, dust grain 
properties may depend on metallicity and on the age of the 
stellar populations. Stellar ages at cosmological distances 
have a strong bearing on the origin of the dust (stellar 
winds of evolved stars versus supernovae). 

In this paper, we present global fluxes from 0.3 /im to 
90 cm for 7 BCDs with metallicities ranging from 1/50 
Zq to solar. In addition to our new near-infrared (NIK) 
images, we have also acquired data from the HST and 
JCMT/SCUBA archives, and gleaned optical, infrared, 
and radio data from the literature. In Sect. |21 we de- 
scribe the sample, the data, and the data reduction. Sect. 
|2|presents the SEDs for the sample BCDs, and Sect.01dis- 
cusses their interpretation in terms of "active" and "pas- 
sive" star formation. New SED models for the BCDs in 
our sample with sufficient data are given in Sect. El and 
scaling laws among mid-, far-infrared, and radio emission 
are described in Sect.|B| We discuss the implications of our 
results on star-formation history and metallicity in Sects. 
ElandlHl and in the context of high-redshift star formation 
in Sect. El If low-metallicity BCD SEDs differ from the 
usual templates (e.g., M82 and Arp220), it will be nec- 
essary to reassess the interpretation of the faint sources 
at z ^ 6 expected to be detected with the new sub-mm, 
mm, and radio facilities such as APEX, ALMA, SKA, and 
LOFAR. 

2. The sample and the data 

Our sample is based on all the low-metallicity BCDs avail- 
able from the JCMT/SCUBA archive as of December, 
2003. This resulted in seven BCDs with metallicities rang- 
ing from solar to ;$ 1/50 Zq. Our sample includes the 
two most metal-poor star-forming galaxies known, I Zw 18 
( < 1/50 Zq) and SBS 0335-052 (1/41 Zq). Tableffllists 
the basic properties of the sample. Properties are taken 
from NED, with the exception of the metallicities, which 



were gleaned from the references listed in the last col- 
umn of the Table. Distances were calculated assuming a 
Hubble constant Ho — 70kms~^ Mpc~^, after correct- 
ing the heliocentric radial velocities to the CMB reference 
frame according to the prescription of the Third Reference 
Catalogue (RC3, de Vaucouleurs et al. 1991). While this 
sample is neither complete nor homogeneously selected, it 
contains many heavily studied BCDs and enables a com- 
parison of their spectral properties. It also has the virtue 
of spanning a factor of 50 in metallicity, and a factor of 
/t 10^ in blue luminosity. With the exception of Mrk33, 
all the sample galaxies are BCDs according to the criteria 
(lu minosity, central color, p eak surface brightness) set out 
bv iCil de Paz et al.l ^Q^. Mrk33, with an absolute B 
magnitude of —18.4, is slightly too luminous to be consid- 
ered a true dwarf galaxy. 

2.1. SCUBA data reduction 

Data reduction of SCUBA data was carried on in the 
standard way, using the dedicated software SURF by 
jjenness fc Lightfootl l)l998j) . The data were first flat- 
fielded to take into account different bolometer sen- 
sitivities, then corrected for atmospheric extinction. 
Atmospheric opacities were mainly derived from the 
225GHz opacity monitor, using the well studied relations 
between the optical depth at 225GHz and those at 450/Ltm 
and 850/j,m. Bolometers that appeared free of source emis- 
sion were used to remove residual sky noise fluctuations. 
Flux calibration was achieved from observations of sec- 
ondary calibrators or Uranus, when available. 

Five objects were observed in jiggle-map mode, with 
data available for both the short (450/im) and long 
(850/im) wavelength arrays. Maps were produced after re- 
binning the data in the RA-Dec plane. Total fluxes were 
measured integrating the source signal (an d the calibra- 
tor) o ver a s elected apertu r e, as d escribed in lDunne et alJ 
(200(f} and Ijenness et alJ l|2002|) . For detected sources, 
apertures were selected to enclose regions where S/N>2 
(unless they were smaller than a 40" circular aperture, 
which was used in these cases) . II Zw 40 was observed in 
10 nights between February 2000 and January 2001, for a 
total number of 315 integrations (10 integrations require 
about 20 minutes to be completed, with the bolometer be- 
ing on-source half of the time) . The sky opacity at 850/im 
was in the range Tgso — 0.2 — 0.3. The object was detected 
both at 450/im and 850/im. At the longer wavelength 
the galaxy is well resolved by SCUBA (FWHM«15"). 
He 2—10 was observed in two nights in December 2000 (95 
integrations) with sky opacity rgso = 0.2 — 0.3. The ob- 
ject has also been detected at both wavelengths. II Zw 70 
and NGC 5253 were observed together in a single night 
(1999, January 13) with rgso — 0.3 — 0.4 (20 integrations 
each). Only NGC 5253 was detected, at 850/im. Finally, 
Mrk 33 was observed during several nights between March 
1998 and January 1999 (106 integrations), with rgso in the 
range 0.1 — 0.5. It also has been detected at 850/im. Fluxes 
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Table 1. Sample galaxies 



Name 


a 


b 


z 


D (Mpc)" 


Mb 


12 + log(0/H) 


Zq/Z" 




Ref.'* 


IZwl8 


0.4 


0.2 


0.0025 


13.0 


-14.1 


7.19 


55 


0.138 


1 


SBS 0335-052 


0.2 


0.2 


0.0135 


55.7 


-16.7 


7.32 


41 


0.202 


1 


IIZw40 


0.6 


0.2 


0.0026 


12.1 


-14.9 


8.09 


7 


3.538 


2,3 


IIZw70=UGC9560 


0.7 


0.3 


0.0041 


19.5 


-16.7 


8.06 


7 


0.053 


3 


NGC 5253 


5.0 


1.9 


0.0014 


4.1'= 


-17.2 


8.19 


5 


0.242 


1,3,4 


Mrk33=Haro2 


1.0 


0.9 


0.0049 


23.4 


-18.4 


8.40 


3 


0.052 


5 


He 2-10 


1.9 


1.4 


0.0029 


9.0-'^ 


-17.3 


8.93 


1 


0.481 


4 



"Distances calculated with Ho = 70kms ^ Mpc ^, and correcting to the CMB reference frame as described in the text. The 

distance for I Zw 18 derived bv lOstlinI fcOOOll is 12.6 Mpc 

'Assuming solar 12 + logO/H = 8.93 jAnders fc Grevess3ll98<ih 
■^Galactic extinction Ab (mag) taken from Schleeel et al,| jl99?i) . 

''R eference for met allicity: (1) Izotov fc Tfauan (1999ll : (21 iGuseva et all J2000l) : (3) iKobulnickv fc SkillmaJ (1199611 : 

f4l lKobulnickv et alJ (il99fll ): f5) lLegrand et alJ l)l997fl . 

"^Distance taken the Cepheid variable estim ate given bv lSandaee ieiandliii). 
■'^Distance taken from IKobulnickv fc Johrison (19991) ■ 



for these objects (and 2-a upper limits) are reported in 
TablesEl To our knowledge, the data for only two of these 
objects have be en previou sly published (namely NGC 5253 
and Mrk 33 .James et alJ [2002^: the fluxes from the litera- 
ture are in excellent agreement with our independent data 
reduction, the results of which are reported in Tabled 

The two remaining objects of our sample, I Zw 18 and 
SBS 0335—052 have been observed in photometry mode 
in February 2000, in nights with optical depths rgso = 
0.3 — 0.4. A total number of 497 integrations were spent 
on IZwl8 and 200 on SBS 0335-052. Neither object was 
detected; upper limits are given in Tables |21 

2.2. Data from the literature 

We searched the literature for optical, NIR, mid- and far- 
infrared, sub-mm, mm, and radio data. Data from our 
own observations were incorporated for Mrk 33 (Haro2), 
IIZw40, and SBS 0335— 052. Because we aim in this pa- 
per to derive the SED for the entire galaxy, rather than 
individual regions, we have placed greater emphasis on 
the largest aperture available and total fluxes. This is a 
problem particularly for the older mid-infrared (MIR) ob- 
servations with photometers, which were, of necessity, per- 
formed with small apertures. Global fluxes are frequently 
unavailable, and comparison of different aperture sizes at 
different wavelengths is the major source of the uncertain- 
ties in our SEDs. 

For wavelengths < 5 /im, Galactic exti nction was cor- 
rected according to the values given by ISchleeel et alJ 
(|l998^ (see Table ^ , a nd using the Galactic extinction 
curve with Ry = 3.1 bv lCardelli et alJ (|l989l) . Corrected 
optical m agnitud e s were converted to fluxes using the zero 
points by'Besselll (|l97i), and NIR/MIR ones given by the 
UKIRT web page^ The UKIRT web page is based on the 



calibr ations given bv lBeckwith et al.l(ll97fil) and lTokunagal 
()l984}) . Interpolation was applied as necessary. 

Fluxes from the Infrared Astronomical Satellite 
(IRAS) were taken either from the NASA/IPAC 
Extragalactic Database (NED), or from published papers. 
Although data were examined on a case-by-case basis, 
single-dish radio observations were generally preferred to 
high-resolution VLA data, because of the possibility of 
losing flux from low spatial frequency filtering. 

2.3. ISO data for SBS 0335- 052 

To re-evaluate the 6 5 um flux for SBS 335—052, in light of 
the Spitzer^ results ( Houck et al.l20o4) . we extracted from 
the ISO archive the ISOPHOT data for this object. The 
pipeline reduction and calibration was used, and aperture 
photometry was performed with IRAF^ . We obtain a weak 
(~ 2a) detection for SBS 0335-052 at 65 /Ltm of ~44mJy. 
This is substantia ll y low er than the 112 mJy estimate by 
IPlante fc Sauvagj |2002l l. but in much bett er agreement 
with t he Spitzer IRS spectrum reported bv iHouck et alJ 
ll20fl4l). 



2.4. NIR data for Mrk 33 and IIZw40 

For two BCDs in our sample, we present new NIR imaging 
data. J (1.2 /im), H {1.6 fim), and K {2.2 fim) images of 
IIZw40 were acquired with the 3.8-m United Kingdom 



^ www . j ach . hawaii . edu/ JACpublic/UKIRT/astronomy/ 
conver .html 



^ The Spitzer Space Telescope is a NASA mis- 
sion managed by the Jet Propulsion Laboratory. See 
www . spitzer . caltech . edu. 

IRAF is the Image Analysis and Reduction Facility 
made available to the astronomical community by the 
National Optical Astronomy Observatory, which is operated 
by AURA, Inc., under contract with the U.S. National Science 
Foundation. 
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Table 2. Total SCUBA fluxes for sample galaxies 



Name 


450/im 


850/xm 


Comments 




(mjy) 


(mJy) 




IZwl8 




< 2.5 




IIZw40 


1550 ±350 


90 ±10 


Resolved at 850//m 


SBS 0335-052 




< 5 




IIZw70=UGC9560 




< 20 




NGC 5253 


<2200 


180 ± 20 


Resolved at 850/im 


He 2-10 


580 ± 150 


140 ±15 


Resolved at 850/im 


Mrk33=Haro2 


<170 


50 ±10 


Resolved at 850/im 



Infrared Telescope (UKIRT^) equipped with IRCAM3. 
The observations are part of our ongoi nK project of NIR 
imaging and spectroscopy of BCDs (see lHunt et aljEoO^) . 
The IRCAM3 plate scale is 0'.'28 per pixel, with a total 
field-of-view (FOV) of 72" x 72". We acquired J, H, and 
K images of Mrk33 with ARNICA mounted on the 1.5-m 
Italian Telescope at Gornergrat (TIRGO^). The ARNICA 
plate scale is 0"97 per pixel with a FOV of 4.1'x4.1'. Both 
cameras house NICM0S3 256^ detectors. 

Each galaxy was observed by alternating frames on the 
source and on adjacent empty sky positions, beginning and 
ending each observing sequence with a sky position. For 
the UKIRT data, before the beginning of each sequence, 
dark exposures were acquired with the same parameters 
as the subsequent science frames. Individual frames were 
dark-subtracted (in the case of IRCAM3) and flat-fielded 
with the average of adjacent empty sky frames, after edit- 
ing them for stars (to avoid "holes" in the reduced frames) 
and applying an average-sigma clipping algorithm. The 
reduced frames were then aligned and averaged. All data 
reduction was carried out in the IRAF environment^. 

Photometric calibration at both telescopes was per- 
formed by obs erving standard stars fro m the UKIRT Faint 
Standard List teawarde n et al. before and after the 

source observations. Each standard star was measured in 
several different positions on the array, and flat-fielded 
by dividing the clipped mean of the remaining frames 
in the sequence. To correct the standard-star photometry 
for atmospheric extinction, we used the mean extinction 
coefficients given by each observatory. Virtual aperture 
photometry was performed for both objects, but Mrk33 
is sufficiently large that we opted for asymptotic magni- 
tudes derived from growth curves, rather than truncating 
to a given aperture. The magnitudes derived in this way 
agree to with i n 0.1 mag to the total magnitudes given by 
iJarrett et al Galactic extinction was corrected for 



* The United Kingdom Infrared Telescope is operated by the 
Joint Astronomy Centre on behalf of the U.K. Particle Physics 
and Astronomy Research Council. 

^ TIRGO (Gornergrat, CH) is operated by the Istituto di 
Radioastronomia-Sezione Firenze. 

® IRAF is the Image Analysis and Reduction Facility 
made available to the astronomical community by the 
National Optical Astronomy Observatory, which is operated 
by AURA, Inc., under contract with the U.S. National Science 
Foundation. 



and conversion to flux was achieved as described in Sect. 

ma 

3. Spectral energy distributions 

The SEDs of our sample galaxies are presented graph- 
ically in Fig. n and Fig. |5| The data are given in the 
Appendix as separate tables, together with the references. 
Although we call them "SEDs" , we have plotted here the 
flux (mJy) versus wavelength, a usual, but conceptually 
unsound, convention. We have done so because it is much 
easier to intuitively interpret the trends in the distribu- 
tions, and identify low flux levels because of the vertical 
scaling common to all graphs. 

We have compared the SEDs shown in Figs. Q and [21 
with the s pectrophotometric model s for M 82 and Arp 22 
given by ISilva et alJ l)l998|) and iBressan et all l)2002|) . 
These two galaxies were chosen because they are fre- 
quently used as t emplates for photom etric redshift de- 
terminations fe.g.. iHughes et aLll2002() . The models have 
been normalized by calculating the mean offset over the 
entire ensemble of data points in logarithmic flux units. 
Such a procedure is equivalent to adjusting the vertical 
scale so as to minimize the residuals data— model. The 
global spectral normalization is crucial to interpret devi- 
ations from the models because of the insidious effects of 
dust extinction, radio free-free absorption, and different 
dust temperatures. 

Problems with total versus large- aperture flux are ap- 
parent in NGC 5253 and He 2-10 (Fig. HI. In both cases, 
excluding the lower-flux points does not change the nor- 
malization relative to the models; we have shown all data 
and included them in the best-fit assessment (see below). 

3.1. Deviations from the standard starburst templates 

To assess how weU the SEDs of M82 and Arp 220 fit 
the low-metallicity BCDs in our sample, we calculated 
the average RMS (logarithmic) residuals after normal- 
ization. We also divided the SED into four spectral re- 
gions, and calculated the average mean offset (normal- 
ization) and RMS residuals individually for each region: 
optical±NIR (A < 4/im), MIR (4 < A < 30 ^im), FIR-mm 
(30 < A < 7000 ^m), and radio (A > 7.0 mm). The models 
were interpolated to the wavelength of the observation. In 
every case, the model which gave the smallest RMS resid- 
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Fig. 1. Spectral flux distributions for He2-10, II Zw 40 NGC 5253. and SBS 0335-052. Upper limits are sliown as 
open symbols. Also plotted are the standard templates l)Silva et alJIlQQsl) for two prototypical starbursts, M82 and 
Arp 220. 



uals after normalization is M82; Arp 220 is not a good 
approximation to any of the galaxies in our sample. 

The results of the normalization/RMS procedure are 
given in TableHfor the best-fit model (M82). Column (2) 
gives the (logarithmic) RMS residuals over the entire SED, 
and Columns (3) through (6) contain first the mean offset, 
and then the RMS residual for the individual spectral re- 
gions. In all columns, the number enclosed in parentheses 
is the number of data points. The normalization procedure 
guarantees that the global mean offset is zero (see previ- 
ous section). All units in the table are logarithmic. Figure 
shows the deviations from the M 82 SED graphically. 



Inspection of Figs. El and Table El reveals several 
important features. 

Radio: Three of the BCDs (IIZw40, NGC 5253, 
SBS 0335—052) have approximately flat global 
radio emission; on galaxy scales the radio emission 
appears to be dominated by thermal processes. 

FIR-mm: In three objects (II Zw 40, NGC 5253, He 2-10), 
the 60-100 i-iTu emission is significantly larger than the 
best-fit template. In SBS 0335— 052, the 65 /im point 
falls well below it. 

MIR: In the objects for which we have considerable 
spectral coverage at the mid-infrared wavelengths 
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Fig. 2. Same as in Fig. H] but for II Zw 70, I Zw 18, and Mrk 33. 



(IIZw40, NGC5253, SBS 0335-052), the observed 
MIR spectrum is not well fit by the M82 template. 
First, the flux for Ax!.12^m exceeds the template 
significantly. Second, there are no PAH features (or 
Aromatic Features in Emission, AFEs) in these ob- 
jects, unlike M82 (and Arp 220). 
Optical+NIR: Three BCDs (IIZw40, IZwl8, 
SBS 0335—052) have optical+NIR emission which 
falls well below the M82 template. IIZw70 appears 
to have an optical excess relative to M82, but a NIR 
deficit. 



4. Active and passive star formation 

The preceding section suggests that three galaxies deviate 
from the M 82 SED over the entire spectral range covered 
by our SEDs: IIZw40, NGC5253, SBS 0335-052. IZwl8 
differs significantly in its optical+NIR emission, but the 
mid- and far-infrared coverage is too sparse to make an 
accurate comparison; this situation should be remedied 
by Spitzer. The remaining galaxies, Mrk 33 and IIZw70, 
follow the template reasonably well. 

The deviant galaxies, IIZw40, NGC5253, and 
SBS 0335— 052, share several characterist ics, namely: 
(i) thev host Super Star Clusters (SSCs) ijThuan et alJ 
Il997t ICalzetti et all ^9?; 'Bec k et all l2002l) : (ii) unlike 
most galaxies, their FIR emission peaks at or shortward 
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Table 3. Deviations from the starburst template, M82 



Name 
(1) 


Global RMS 
(2) 


Optical+NIR 
(3) 




MIR 

(4) 




FIR-mm 

(5) 






Radio 
(6) 




IZwlS 


0.44 (20) 


-0.11 


0.57 


(8) 














-0.01 


0.24 


(10) 


IIZw40 


0.66 (47) 


0.31 


0.79 


(20) 


-0.37 


0.48 


(13) 


0.06 


0.31 


(6) 


-0.20 


0.28 


(8) 


SBS 0335-052 


0.43 (106) 


-0.05 


0.62 


(8) 


0.14 


0.37 


(86) 


-0.36 


0.36 


(3) 


-0.13 


0.35 


(6) 


II Zw70=UGC9560 


0.34 (25) 


0.19 


0.38 


(10) 


-0.37 




(1) 


-0.16 


0.19 


(2) 


-0.15 


0.24 


(11) 


NGC 5253 


0.35 (48) 


0.39 


0.46 


(12) 


0.14 


0.20 


(23) 


0.12 


0.13 


(4) 


-0.19 


0.26 


(8) 


He 2-10 


0.29 (41) 


0.20 


0.29 


(16) 


-0.05 


0.12 


(6) 


-0.09 


0.37 


(5) 


-0.17 


0.18 


(14) 


Mrk33=Haro2 


0.25 (29) 


0.13 


0.24 


(13) 


-0.19 


0.08 


(2) 


0.02 


0.15 


(3) 


-0.11 


0.24 


(10) 



■a 

K 



-1 




«• NGC5253 

O IIZw70 

• nZw40 
□ IZwlS 

■ SBS0335-052 

* Mrk33 
A Hc2 10 

I , I 

Optical MIR 



I , I 

FIR-mm Radio 



Fig. 3. The mean logarithmic residuals of each galaxy rel- 
ative to the best-fit starburst template, M 82 as a function 
of spectral region. The exact horizontal position is arbi- 
trary, as they have been slightly offset for better visibility. 
Error bars are the error in the mean, so divided by ^/n 
relative to the values in Table |31 



of eO/xm l|Vader et all IPlante & Saiivagel 120021) : 

(m) there are no PAH features in thei r MIR spec- 
tra llRoche et al.lll99ll: iThuan et al]ll999bt iHouck et alJ 
20041 ): (iv) theyhave K-L > 1 llMoorwood fc Glassll982t 



Jov fc Leste'^ligSStlHunt et alJl200l|) : and ( y) their radio 



spectra are risiiig on small spatial scale s ([Turner et alJ 
l200nHReck erani2002HHunt et a,lJl2004ah . The first four 
characteristics are associated with what has been called an 
"active" mode of star formation (e.g.. iHunt et al.ll2002l: 
iHirashita fc HuniJl2004l) . BCDs which form stars in this 
way have star-forming regions which are small ( ^ 100 pc) 
and dense (ug >, 500 cm~^). Star formation which oc- 
curs in larger, more diffuse regions was dubbed "passive" . 
Such a "dichotomy" does not occur only in low-mass low- 
luminosity BCDs; there is evidence that star formation is 
"bimodal" even in giant infrared-luminous galaxies, with 
more intense starbursts tending to be more heavily ob- 
scured and compact (,Takagi et al...2003i|"l . 



"Active" regions tend to be more efficient at shielding 
ultraviolet (UV) radiation because of the greater concen- 
tration of dust. Hence the optical depth of the UV pho- 
tons becomes large, and dust reaches high er temperatures 
than it would in more diffuse regions ijHirashita et alJ 
Because of the high stellar concentration and 
low metallicities in su ch BCDs, the UV radiation field is 
also very intense (e.g.. lVerma et al.ll2003), wh ich tends to 
suppress PAH emission l|Puget &: Leged 1989|) . "Passive" 
star-forming episodes do not engender SSCs, although 
they can generate massive stellar clusters; dust tempera- 
tures are lower, and the MIR spectrum shows conspicuous 
PAH features. The deviant galaxies, IIZw40, NGC 5253, 
SBS 0335-052, can all be classified as "active" BCDs; they 
are also characterized by compact dense star-forming com- 
plexes, as judged from HST images and optical spectra 
(Hunt & Hirashita 2004, in preparation). 

It is likely that the fifth trait shared by the deviant 
galaxies -compact radio sources with rising spectra- is also 
associated with active star formation. The theoretical ba- 
sis for the active/passive "dichotomy" lies in the size and 
densi ty of the star-forming complex ijHirashi ta & Hunt! 
l2004l) : compact dense regions form stars "actively" and 
larger more diffuse ones "passively" . Rising (inverted) 
radio spectra surrounding individual stars, indicative of 
dense o ptically thick Hii regions , are abundant in the 
Galaxy llWood fc Churchwelll989l) . but in external galax- 
ies have been observed in associa tion only with compact 
( ;$ 10 pc) massive star clusters llKobulnickv fc Johnsor 
1999[ iBeck et all l2000t iTurner et alJ l2000t iBeck et al. 
2002|) . The inferred densities in these small regions are 
extremely high, typically ^ 2000 cm"'^. Hence, the two 
properties required for a star-forming complex to be clas- 
sified as "active" are present in these inverted spectrum 
sources, and they may be an integral part of the phe- 
nomenon. 

In addition to the three obvious "deviants" , there may 
be another active BCD in our sam ple: He 2—10. It contains 
both SSCs (^Johnson et al."2000')_^ an d risin g spectrum ra- 
dio sources (^Kobulnick v fc JohnsonI [19991) . However, its 
FIR dust temperature is similar to normal spiral s, and its 
K ~ L color is ~ 0.8 ("Gla ss fc Moorwoodlll985l). smaller 
than the typical values of "active" starbursts ()Hunt et alJ 
|2002|) . Its MIR spectrum also contains PAH features, to- 
gether with a deep 10 /im absorption trough ((R,oche et alJ 
il991i) . Hence, He 2— 10 may be active on small spatial 
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scales, but globally "passive". At solar metallicity, this 
galaxy is the most metal-rich in our sample. Its relatively 
high metal abundance may imply a more evolved ISM 
which would tend to dilute "active" features and lower 
the dust temperatures observed over large regions. 

It is as yet unclear why some regions form stars 
actively and others passively. Metallicity appears to 
have little or no influence sinc e IZwl8 forms stars 
passively l)Hirashita fc Hunt! l2004l) . while SBS 0335-052 
forms them actively, d espite a similar (low) metal abun- 
dance. iHirashita fc Hu nt ( 200^ speculate that high pres- 
sure and shock compression may be affect the bifurca- 
tion into active and passive regimes. Star-formation rate 
(SFR) may also be important, since active regions tend to 
have higher SFRs than passive ones; however this could 
be an effect rather than a cause. Active regions are also 
invariably young ( ^ lOMyr): indeed optically thick Hii 
regions as signified by warm dust emission and rising ra- 
dio spectra app ear to be the first stage of intense star 
forma tion (e.s.. iBeck et al.lll99fil: iKobulnickv fc .TohnsoiJ 
ll999HVacca et al.l2002|) . For a region to be defined as "ac- 
tive" , it may be necessary to observe it in its early stages. 
Passive star formation on the other hand can be either 
young or more evolved, and both regimes may be associ- 
ated with superbubbles or expanding shells. Indeed, with 
the exc eption of IIZw7 which is a probable polar-ring 
galaxy ijCox et al.ll200lh. all of the galaxies in our sample 
contaiii such s t ructures llMartirJl99a:lMarlowe et al .Il997t 
Martini Il997t iMendez fc EstebanI Il997t ISummers et all 
200l|) . Also with the exceptio n of II Zw70, the sample 
galax ies host Wolf-Rayet stars llCnntilll 99ll |Cuseva et a,]" * 



l200nl) . implying that the star-formation episodes must be 
younger than a few Myr ijSchaerer fc Vaccalll998|) . 



4.1. Optica l/N I R colors 

The hot dust that is associated with active star forma- 
tion is revealed in t he K — L color, an d is usually already 
evident in H - K llHunt et al.ll20ol ^. However, optical 
and NIR colors by themselves may be an ambiguous diag- 
nostic of the active/passive distinction. This point can be 
clarified by examining the colors of our sample galaxies. 

The BCDs which we have classified as "ac- 
tive" (NGC5253, HZw,40, SBS 0335-052, and possibly 
He 2—10) are characterized by a wide range of optical/NIR 
colors. The "passive" galaxies on the other hand (Mrk33, 
IIZw70, IZwl8) have uniformly blue to extremely blue 
optical/NIR colors. This is shown in Fig. 0] where opti- 
cal/NIR colors are plotted for the sample galaxies. Also 
shown in the figur e are the Starburst99 models for an in- 
stantaneous burst l|Leitherer et al.ll999t SB99), with 1/20 
Zq to Zq] the SB99 predictions include nebular contin- 
uum emi ssion. The grid gives the observed range for spiral 
galaxies ijde Jondll996|) . The Zq curve is the one which 
touches the spiral galaxy grid at roughly 10 Myr because 



of the onset of the asymptotic giant branch. Fig. 0| shows 
that IZwlS, SBS 0335-052, IIZwTO, and Mrk33 are sig- 
nificantly bluer than normal spirals. Only the solar metal- 
licity SB99 models are roughly consistent with the col- 
ors (although none o f them predicts the colors of I Zw 18, 
see iHunt et aLlEool . and these only for the ages of the 
models, 1 Gyr or less. The remaining galaxies, NGC 5253, 
IIZw40, and He 2— 10, are much redder; they are only 
slightly bluer (NGC 5253) or consistent with the colors of 
early spirals. Extinction, the effect of which is shown as 
an arrow in the Fig. 0] would bring the colors blueward, 
similar to those of the "passive" BCDs. Therefore, pas- 
sive BCDs seem to be associated with significantly blue 
colors from 0.4 to 2 /im, but the colors of active BCDs 
can vary substantially. Part of this variation may be due 
to dust reddening; part of it may be an age effect. With 
optical/NIR colors alone, this ambiguity is difficult to re- 
solve. 

The bluish or normal colors of the active BCDs are 
in apparent contradiction with the substantial emission 
at A > 2 fim on a global scale; this is particularly 
striking in the case of SBS 0335— 052, the prototype of 
an active BCD. The seeming paradox can be under- 
stood through optical depth effects and the morphol- 
ogy at different wavelengths. The dust extinction is so 
high on small scal e s in SBS 0335— 052 (Ay ^ 12 mag, 
iThuan et all Il999bl IPlante fc Sauvagell2002() that more 
than half the star formation is hidden at 2 /xm l|Hunt et alJ 
l200l|) . Thus, in the optical, we are seeing only the outer 
shell of the star-forming regions, the inner parts of which 
are hidden by dust. How much of the r egion is visible de- 
pend s on the geometry of the dust (e.g.. lPlante fc Sauvagj 
|2002|) . and on the configuration of the emitting regions. In 
the optical, extended ionized gas contributes significantly 
to th e global emission in SBS 0335— 052 ijTh^ian et alJ 
I1997I) . while in the infrared, the emission is extremely 
compact, confined to the most luminous star clusters 
l|Dale et al.ll2001btl . 

A mixture of these effects is probably responsible for 
the wide range of colors in what we have called "ac- 
tive" BCDs. In NGC 5253, a single cluster dominates the 
SED beyond 3 /im, while the entire galaxy contributes to 
the optical and NIR emission (jVanzi fc Sauvagc 2O0I . 
The same is essentially true for IIZw40 and He 2— 10 in 
which a large fraction of mid-infrared and radio emis- 
sion comes from compact star clusters in a region less 
than 200 PC in diameter llKobulnickv fc .TohnsonI Il999l: 
' 02). 



However, if not accompanied by a bluish J — H, then red 
H — K IS probably better ascribed to extinction effects. 



iBeck et alJ2'00ltlVacca et a l.'2002':' Beck et al .'200?). Like 
SBS 0335—052, in the se objects there is substantial extinc - 
tion {Av ^ 10 mag, iKawara et al]ll989HHo et alJll990|) . 

Dust extinction is generally much lower in diffuse extended 
regions than in com pact heavily obscured star clusters 
llCalzetti et al.lll997|) . 

5. Modelling the dust in active BCDs 

The notion that low metallicities imply low dust content in 
BCDs is apparently dispelled by our data. If anything, the 
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B-R 

Fig. 4. Optical and optical/NIR colors of the sample 
galaxies. B — K versus V — R are shown in the upper 
panel, and B — H vs. B — R'm the lower. I Zw 18 is shown 
as an open (blue) square; II Zw 70 as an open (blue) circle; 
Mrk33 as an (blue) asterix; He 2— 10 as a filled (green) 
triangle; SBS 0335-052 as a filled (red) square; IIZw40 
as a filled (red) circle; NGC 5253 as an open (red) six- 
point star. The insta ntaneous-burst Starburst99 models 
llLeitherer et al for Zq/ 20 to Zq are shown as (ma- 

genta) solid lines (from 1 Myr to 1 Gyr), a nd the observe d 
colors of spiral galaxies as a (cyan) grid l)de Jondll996^ . 
The arrow in both panels shows Ay = 1 mag, acco r ding t o 
the interstellar extinction curve bv lCardein et all lll989t) . 



active BCDs in our sample have a larger dust fraction rel- 
ative to the radio continuum than M82. This is true both 
for warm dust at MIR wavelengths, and for cooler dust in 
the FIR/mm. Our new SCUBA data allow a more com- 
plete estimate of the dust mass in these BCDs, but first we 
need to fit the SED. To do this we have adopted DUSTY, 
a code which solves the radiation transport problem in 
a dusty environment ((Ivezic fc Elitzur..l997; Jvezic et alj 



The solution is based on a self-consistent equa- 
tion for the radiative energy density, and includes the ef- 
fects of dust scattering, absorption, and emission. Such 
a treatment is important in the active BCDs we study 
here because of optical depth effects; simple screen models 
do not properly account for dust self-absorption. Because 
the code assumes thermal equilibrium, stochastic emis- 
sion from very sma ll grains (VSGs) is not considered (e.g., 
iDesert et a,1.lll 99n|) . While this is not a significant problem 
for dust mass estimates because of the small contribution 
of such grains to the total dust mass, it may affect the 
parameters we infer from the SED fit such as temperature 
and dust distribution. 

Because the active BCDs in our sample all host com- 
pact SSCs or ultradense HII regions or both, we have as- 
sumed a spherical shell geometry for the dust emission. 
The inside of the shell is exposed to radiation from a young 
stellar population at the appr opriate metallicity and age. 
In general, for any geometry l|lvezic fc Elitzuilll997l) . the 
scaling properties of the radiative transfer problem leads 
to SED families which, for a given grain type, are governed 
by the spatial density distribution of the dust; the optical 
depth r deter mines the position within the family. As first 
pointed out bv iRowan-RobinsonI l)l980|) . neither the over- 
all luminosity nor absolute scales of density or distance are 
required for the exact solution; hence, the problem can be 
defined with relatively few quantities. With DUSTY, it is 
necessary to specify the temperature Tin at the innermost 
boundary of the shell, and the thickness of the shell in 
terms of the ratio of the outer to the inner radius, i^ut (= 
Rout/Rin)- For simplicity, we have used one-zone models, 
so it is necessary to specify only a single power law index 
p for the radial run of the mass density, R~p. 

DUSTY currently only supports single- type grains, 
namely a single size and chemical composition. However, 
it allows the simulation of more realistic grain populations 
through the appropriate weighted means. Thus, to define a 
grain population, .Dt/STy requires the minimum and max- 
imum grain size, amin and flmax, assumed to be the same 
for all chemical species; the power law index q describing 
their size distribution, a^"^; and their chemical compo- 
sition. Several species are supported by D USTY includ- 
ing c lassical silicate and graphite grains llDraine fc Lee 
19841) and amorphous carbons (AMCs) (e.g.. iHannei 
1988t IZubko et al.lll99'^ . The fractional number abun- 



dances must be specified for each grain type. 

We implemented the SED fitting of our active BCDs 
by performing a grid search over various dust density dis- 
tributions, l^ut, dust compositions, and Tin. This last im- 
poses the physical size of the innermost shell boundary 
i?i„ because of its dependence on the intensity of the cen- 
tral source and the normalization of the radial coordi- 
nate. We defined three basic dust compositions for th e 
grid search: (i) a "standard" one ijDraine Sz Leel[l98^ : 
(ii) a mixture of silicates, graphite grains, and amorphous 
carbons fe.g.. lWeingartner fc Drainel200lf) : [iii) a silicate- 
free mixture with only graphite grains and amorphous car- 
bons. Grain sizes were assumed to be distributed accord- 
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ing to a power law, with the index q f ixed t o —3.6, in 
accordance w ith IWeingartner fc Draind ll200ll). an d very 
similar to the lMathis. Rumpl. fc Nordsieckl l)l977l MRN) 
value of —3.5. Grain sizes were taken to be "standard" 
MRN, ranging from 0.005 /xm to 0.25 /im. The central ex- 
citing source spectrum was taken from SB99 for the ap- 
propriate age {'-^ 3 Myr, see Sect. 15. 1|) and metallicity as 
determined from the optical spectra. For each metallicity 
(1/20 and 1/5 Zq), more than 2000 DUSTY uYodeis were 
generated in the grid, with power-law indices p ranging 
from to 1, y^ut from 20 to 1000 logarithmically, and Tin 
from 400 to TOOK, tv was varied logarithmically from 0.1 
to 100. 

Because of the importance of thermal radio emission 
to the observed SED, we first fit a single power law to 
the radio spectrum for A > 1 cm. This fit was subtracted 
from the observed spectrum for A < 1 cm to obtain the 
pure dust contribution. The best-fit DUSTY model was 
determined by minimizing the logarithmic residuals of the 
radio-subtracted dust spectrum and the _D[/;ST'y models 
from 5 /im to 1cm. Once the best-fitting DUSTY model 
was achieved in this "coarse" grid, the parameter space 
of the models was refined to describe the SED as well as 
possible. This resulted in typically another 20—50 mod- 
els being generated. With the best-fit parameters, we also 
investigated the effects of varying the age of the stellar 
population and the grain size distribution and chemical 
composition. 

5.1. DUSTY model results 

The best-fitting £'[/S'TF models are shown in Fig. |S1 and 
the radial run of temperatures in Fig. (H] The dust pa- 
rameters and radial distributions are reported in Tabled] 
Uncertainties in the model parameters are inferred from 
the grid spacing with which we searched the parameter 
space, together with the consequent variation in the mean 
logarithmic residuals. The uncertainties in the tempera- 
ture at the inner shell boundary Tin is roughly 25° K, and 
that at the outer boundary Tout on the order of a few (2- 
5°K). As for the dust mass distribution, the power-law 
index has an uncertainty of a few tenths, and the abso- 
lute values of the radii between 10 and 20%. The dust 
optical depth tv is precise to ^--^30% for the larger values 
(ry = 30) and about 20% for the smaller ones [tv ^ 20). 

The optical and NIR data points were not included 
in the fits because of unknown extinction and geometry 
effects. The stellar populations reflected by these data in 
the global SEDs are almost certainly not subject to the 
high extinction experienced by the central star clusters 
responsible for heating the dust (e.g.,|^alzctti ct al. 1997). 
Indeed, Fig.Oshows that if they were, they would not be 
visible at all at optical wavelengths. 

Also given in Table 0] are the inferred dust masses 
Afdust and infrared luminosity Lir from 5 fim to 1 cm; this 
last is approximately the bolomctric luminosity because of 
dust reprocessing of the incident UV radiation. Tir is eas- 




R (pc) 



Fig. 6. Radial run of dust temperature Tdust as predicted 
by the best-fit T'C/S'Tr models. SBS 0335-052 is shown 
by a solid (red) line, He2— 10 by a long-dashed (green) 
one, IIZw40 by a short-dashed (blue) one, and NGC 5253 
by a dotted (red) one. 

ily obtained by integrating the model over the fitted wave- 
length range, and applying the normalization factor used 
to rescale the £)[/5ry models to the observed SEDs. The 
dust mass is less straightforward to derive from the mod- 
els. Althoug h Mdiist is proportional to the optical depth 
TV (see e.g.. lPlante fc Sauvagel l2002V the former depends 
on the amount of dust and its geometry, while the latter 
is also a function of dust composition, its optical prop- 
erties, and the grain size distribution. Hence, converting 
the optical depths given by the models into dust mass in- 
volves calculating weighted averages of these quantities in 
the same way as in the DUSTY code. We have done this 
as well as possible, but the approximations for the optical 
quantity weighting could introduce an additional 10-20% 
uncertainty in the dust masses. The equations which re- 
late Tv to Mdust are given in Appendix^ The main un- 
certainty in Afdust however is the dust geometry, that is 
to say the radial distribution and truncation radius i?out- 
From the difference between the best fit and the next best 
one, uncertainties in the dust mass are usually 20-30%, 
but can be as high (in the case of SBS 0335—052) as a 
factor of 3. 

The most important result of the fitting procedure is 
the good quality of the fits: with the exception of II Zw 40, 
the mean (logarithmic) residuals are 0.1 or better (see 
Table 0J. This implies that both the assumed shell geome- 
try and the approximation of thermal equilibrium for the 
dust emission may not be unrealistic hypotheses. On the 
other hand, there could be a strong degeneracy between 
the temperature distribution of the dust derived from the 
assumption of thermal emission and the stochastic emis- 
sion from VSGs. Either way the dust masses are likely to 
be unaffected because of the neglible contribution from 
the small-grain population. In IIZw40, the quality of the 
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Fig. 5. Best-fitting DUSTY SED models for the "active" galaxies (upper panels shown as solid lines), together with 
the logarithmic residuals of the fits (lower panels). The SEDs for Arp220 and M82 are also shown, repeated from 
Figures ^ and 13 



fit (mean residuals 0.15) is almost certainly affected by 
aperture effects since the cool-dust emission at 850 /im ap- 
pears extended beyond 4 0". The 230 GHz point of 5.9 mJy 
ijKobulnickv et aDll999l) for He 2—10 is impossible to rec- 
oncile with any fit and was disregarded both in the nor- 
malization and in the calculation of the residuals. 

In general, with the exception of SBS 0335— 052, we 
find that the dust mass distribution is not strongly peaked; 
either uniform or low-index power laws provide the best 
fit for the radial distribution of the dust. Again with the 
exception of SBS 0335— 052, the star clusters modelled by 
our fits have rather thick shells with Rout/Rin^ 500—1000. 
They ar e also rather extended, having radii of ^ 200 pc 
(see also lVanzi fc SauvagQi2004i|') . However, this size is de- 



ceptive, being the truncation radius of the dust shell; the 
dust at these radii is cool, around 20° K (see Table 0] and 
Fig. ini) , and would not be detected in the MIR. The form 
of the SEDs at longer IR wavelengths dictates the presence 
of such cool dust, because of their broad peak at ~60 /xm. 
On the othe r hand, as shown by the IRS spectrum of 
iHouck et al.l l)2004 . the dust in SBS 0335-052 is very 
warm, peaking at ~30/im. It is also very compact; our 
fit gives a diameter of ~ 10 pc, with a dust temperature of 
~ 70° K at this boundary. Our new 65 /im ISOPHOT flux 
of 44mJy is consistent with this, but combined with the 
Spitzer IRS spectrum results in a much different DUS TY 
fit than obtained previously l|Plante & Sauvag3l2002t see 
below) . 
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Table 4. Z^TOTy best-fit parameters for "active" BCDs 



Name 


a 




Tin 


Tout 


Rin 


Ro\it 


TV 


Dust" 




Mdust 






(radio) 




(K) 


(K) 


(pc) 


(pc) 








(Mo) 


(log) 


IIZw40 


-0.5 




600 


18 


0.3 


297 


20 


WD 


3.0 X 10"* 


1.8 X lO'' 


0.15 (19) 


SBS 0335-052 


-0.3 


R-^ 


400 


69 


0.5 


5.4 


30 


WD+Zb 


1.4 X 10^ 


1.1 X 10^ 


0.08 (88) 


NGC 5253 


-0.2 


^-0.3 


450 


19 


0.4 


213 


10 


Grf+Zb 


1.8 X 10^ 


1.6 X 10^ 


0.10 (25) 


He 2-10 


-0.6 




600 


23 


0.4 


200 


20 


WD 


5.6 X 10^ 


9.6 X 10* 


0.09 (10) 



Dust is coded by: WD-|-Hn (0.3 silicates, 0.4 graphite, 0.3 Hanner AMCs); WD-|-Zb (0.3 silicates, 0.4 graphite, 0.3 Zubko et 
al. AMCs); Grf-I-Zb (0.70 graphite, 0.3 Zubko et al. AMCs); WD (0.6 silicates, 0.4 graphite). All grains follow a power-law size 
distribution with q — 3.6 and amin = 0.005 /im, amax ~ 0.25 nm, except for Grf-f-Zb which has amin ~ 0.03 and Omax = 0.3 fim. 
^ Mean residuals from 5 /im to 1 cm. 



The dust in these BCDs also appears to be rela- 
tively "standard" . The same (~MRN) grain size distribu- 
tion and standard composition (0.3 silicates, 0.4 graphite 
grains, 0.3 AMCs), fits all galaxies except for NGC 5253, in 
which there seems to be a paucity o f small grains and sil- 
icates, in agreement with lVanzi fc Sauvage 1,20041 . AMC 
grains are necessary in SBS 0335—052 and NGC 5253, but 
not in He 2-10 and IIZw40 where a "standard" MRN 
composition provides the best fit. 

We found very little if any dependence of the models 
on either the metallicity or the age of the exciting stellar 
cluster. This is probably because luminosity enters into 
the radiative transfer problem only through the flux at 
the inner radius of the shell; once this is specified by the 
temperature Tin and the radius i?in, the problem merely is 
rescaled through normalization. In all cases, we therefore 
adopted a 3Myr cluster for the final best fit, consistent 
with the optical spectra and the presence of Wolf-Rayet 
features (see Sect. 0J), although we used the SB99 stel- 
lar population with the metallicity closest to that of the 
particular BCD. 

5.1.1. IIZw40 

Our model gives a total IR luminosity Lir of 3.0 x 10^ L©, 
in excellent agreement with the OB stel lar luminosity ioB 
of 3 x 10^ Lq rep orted bv iBeck et all (.2002). and with 
earlier estimates l|Wvnn- Williams fc Becklinlll986(l . The 
structure in the starburst of II Zw 40 is however not well 
approximated by a single dusty shell, since high-resolution 
radio maps show four c ompact nebulae separated by 10- 
12 pc l|Beck et al.ll2002|l . The inner radius of 0.3 pc given 
by our best-fit model is probably the best guess for the 
combined emission from these clusters; there is clearly 
warm dust (our model gives Tin — 600 ± 25 K) heated by 
the massive stars. 

Assuming an age of 3-5 Myr, we can use the SB99 mod- 
els to derive the stellar mass responsible for the IR lumi- 
nosity 3. Ox 10^ Lq gives a total stellar mass of 2. Ox 10^ Mq 
associated with the starburst. This is aga i n in v ery good 
agreement with the estimate of lBeck et al.l l|2002t) who find 
a total of 2 X 10® Mq in young stars, and is ^ 10% of the 
total stellar mass estimated bv lVanzi et al.l lll996|) . 

The cool-dust emission is clearly extended in this 
galaxy, since the flux at 40" is about half of the global 



emission in the SCUBA image. Narrow-band images and 
long-slit spectroscopy show that the emission-line region 

in II Zw 40 is more extended than 40" (Vanzi et al 1996|) , 

so the presence of spatially extended cool dust in the 
SCUBA image is not surprising. The diameter of '--^BOOpc 
(where T ~ 20 K) given by our model is much smaller 
than this ^ 10"; this should be however roughly the max- 
imum extent of the dense star clusters that comprise the 
starburst. the cor e alone of which has an extent of ~150 pc 
l(Beck et al.ll2002l) . 

The extinction given by our model of ry = 20 mag is 
significantly higher than previous estimates {tv ~ 3 — 5 ; 
I.Tov fc LesteJl988HWalsh fc R,ovll99.'^IVanzi et alJll99fil) . 
Part of the reason for this may be the assumption of a fore- 
ground screen by these authors. However, the discrepancy 
may also be ascribed to substantial optical thickness; in 
very dusty star clusters, even the NIR hydrogen recom- 
bination lines do not prob e the central most obscured re- 
gions (e.g., SBS 0335^052: iHunt et af]l200lh . 

The best-fit radio index is a = —0.2 which gives large 
£'[/S'TF residuals (cr) — 0.2. We therefore steepened the 
radio spectral index a gradually to —0.5; the residuals 
monotonically decrease to 0.15 as a steepens. We therefore 
adopted a = —0.5 which clearly does not well approximate 
the radio emission at lower frequencies; there is evidently 
substantial free- free absorption in II Zw40 (see lBeck et alJ 

5.1.2. SBS 0335-052 

Our DUSTY ?A for SBS 335-052 is substantially differ- 
ent than that presented bv'Pla nte fc Sauvagd 1*200^ : most 
importantly, it predicts a much lower dust mass and a 
much smaller shell. Our dust ma ss estimate of 1 .1 x 10 ^ Mq 
is a factor of 10 lower than by iHo^ick et al ll2004^■ and 
10'^ tinies lower than previous estimates ( Thuan et alJ 
Il999hl: iHunt et alJl200lt IPlante fc Sauvagell2nn2^ . There 
is roughly a factor of 3 uncertainty on Mdust! our next 
best fit gives a dust mass of ^ 3.7 x 10^ Af© with a 
truncation radius i?out of 11 pc. However, if i?out (hence 
-^^dust) are much larger, the residuals become unaccept- 
ably high. The main reason for the difference is the input 
SED in w hich we have i nclude d the Spitzer IRS spectrum 
given by 'Ho uck et all ([20041) . and our new ISOPHOT 
65 /im photometry (44mJy vs. 112 mJy). As pointed out 
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bv lHouck et al.1 1)2004|) . the spectral turnover at ~ 30 
in the IRS data imphes that the cool dust is warmer than 
previously thought, lowering the inferred dust mass. The 
IR luminosity 1.4 x 10^ L^T) of our model i s also ~3 times 
lower than that of lPlante fc Sauva'g^ l)2002j) , again because 
of the different SED. The definitive answer of how much 
dust is contained in the SSCs of SBS 0335—052 must wait 
for the long-wavelength Spitzer MIPS photometry. 

The other parameters of our best- fit DUSTY model 
are al so rather different from those bv lPlante fc Sauvaeel 
ll2002t) . Although the radial fall-off of their three-zone 
model is the same as our one-zone {R~^) and the optical 
depths are the same (ry = 30), they find a much thicker 
and hotter shell than we do (Tin = 700 K, Yout^lOOO). 
Both of these differences can also be ascribed to the new 
IRS data. The smaller thickness of our model is due to 
the lack of cool dust emission beyond 30 ^m; thick shells 
emit at longer wavelengths because of the greater abun- 
dance of cooler dust. The lower temperature of 400 K is 
dictated by the slope and relative fiux of the IRS spec- 
trum for A £ lO^m; hotter dust (than ~500K) gives too 
much flux, and a slope that is not sufficiently steep. The 
dust parameters we find for SBS 033 5—052 also differ from 
those found by IPlante fc Sauvasrel |2002l) . Although the 
dust composition of both fits is similar, they find evidence 
for a ffatter exponent in the size distribution and larger 
grains, while we find that it is possible to fit the SED with 
"standard" dust sizes. 

Assuming an age of 3-5 Myr, as for IIZw40, we can 
derive the total stellar mass responsible for the IR lumi- 
nosity 2.8 X 10^ Lq and fi nd ^ 10^ A fm, only sfightly lower 
than that found bv .Plante fc Sauvaee 1.2002 ) ■ 2 x 10*^ Mq. 
It i s however much lo wer than the 6.6 x 10^ Mq estimate 
by iHunt et alJ l|200ll) on the basis of Bra. Radio mea- 
surements have shown that the number of ionizing pho- 
tons i nferred from Bra i n SBS 0335— 052 is inconsistently 
high ^ Hunt et a,].l l 2004ah because of the infiuence of stel- 
lar winds; the stellar mass given by yHunt et a l. (200y 
is consequently overestimated. Indeed, embedded massive 
star clusters and SSCs all tend to have masses around 
~ 10^ MfTi LGallaghcr fc Smith 1999: Sm ith fc Gallaghe^ 
l200lt IVacca et al.ri2002tlde Griis et al.ri2004iti r and this 
value is typical of the BCDs in this paper. 

5.1.3. He2-10 

The total Ljr of our best-fit model is 5.6 x 10^ Lq, 
in good agreement w ith the value of 6 x 10^ Lq given 
bv Ivacca et al.1 l|2002fl . Our dust mass estimate Mdust 
10^ Mq is substantially lower than theirs, although con- 
sistent with that of Baas ct al. ( 199^. The stellar mass in- 
ferred from SB99 assuming an age of 3 Myr is 4 x 10^ Mtn, 
consist ent with previou s estimates (e.g., Ijohnson et alJ 
[20 00; Vacca et al.ll2002|) 

IVacca et al.l l)2002l) have modelled the IR emission of 
He 2— 10 using optically-thin shell models of the dust with 
a uniform density distribution. Despite the differences in 



the model assumptions, i?out of our DUSTY model is sim- 
ilar to theirs (100-200 pc for the various emission knots). 
However, the inner and outer temperatures are rather dif- 
ferent: we find Tin = 400 K and Tout = 23 K, while they 
find a factor of two cooler dust at the inner shell bound- 
ary and warmer dust (~ 30 K) at the outer boundary. 
They also find much lower extinction {Ay ^ 8) while we 
find Ty = 2 a va lue c onsistent with the estima te by 
IPhillips et all l|l984l) and iKobulnickv et alJ lll995l). The 
differences in our results and those of I Vacca et al l (|2002h 
are almost certainly due to their assumption of an opti- 
cally thin shell, and a single grain radius of 0.1 fim. The 
high value of Ty given by our mod el is consistent with the 
observation of IVacca et al.1 l|2002fl that "none of the SSCs 
seen in the HST UV and optical images, or the ground- 
based NIR images, is detected in the Gemini iV-band im- 
age, and vice versa" . 

5.1.4. NGC5253 

The IR luminosity of 1.8 x 10^ Lq given by our DUSTY 
m odel is in excel l ent ag reement with that of 1.7 x 10^ Lp, 
of lGoriian eT^ l)200l|i . and of lVanzi fc Sauvagd l|2004li 
after correcting for the different adopted distance. This 
translates into a stellar mass of 1.2 x 10^ Mq, assuming 
an age of 3-5 Myr and using the SB99 models as before. As 
mentioned previously, this age is appropriate for galaxies 
such as NGC 5 253 with Wolf-Rayet features in their op- 
tical spectrum (jSchaerer et al ] ll997l). Our estimate of the 
stellar mass agrees well with that of lCalzetti et al 
and again is typical of the mass of embedded clusters and 
SSCs (see above). 

The SE D of NGC 5253 has als o been modelled with 
L)[/gTFbv lVanzi fc Sauw3 l)2004j) . The dust mass Mdust 
of 1.6 X 10^ Mq given by our best fit is consistent with, al- 
though lower than, theirs (2.3 x 10^ Mq with our Cepheid 
distance of 4.1 Mpc). The geometry of the two fits is also 
similar, with our one-zone best-fit model compared 
to their two-zone distribution, i?""-^. The outer shell 
radii also agree very well, once the different adopted dis- 
tance has been accounted for (216 pc vs. our value of 
213 pc). Moreover, the two dust compositions are simi- 
lar, with a lack of silicates, and 70% (us) to 80% (them) 
in graphite grains. The optical depth ry 10 of our 
best-fit model is in good agreenient w ith that of > 9 
mag measured bv ICalzetti et al.l l(l997^ from optical hy- 
drogen recombination lines, alth ough slightly higher than 
{Av = 8) derived bv lVanzi fc Sa uvagc (2004). 

However, IVanzi fc Sauvagd l)2004|) find substantially 
hotter dust (570 K vs. 450 K) and slightly larger grain 
sizes. The differences b etween our results and those of 
IVanzi fc Saiivagj l(2004l) are almost certainly attributable 
to the differences in fitted grain properties. Our second- 
best fit (also characterized by (cr) = 0.10) is obtained 
with a "standard" grain size distribution (see Sect. |5l and 
Table 0J, and results in a dust mass Mdust of 4.8xlO'^M0, 
with Ty = 8 and Tin = 400 K. This dust mass is more than 
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three times smaller than that obtained with our best-fit 
model which has amin — 0.03/^m (rather than the "stan- 
dard" ainin — 0.005^m). Indeed, the mass of dust with 
given flux from the rmalized grains de pends linearly on the 
grain size (see e.g., lHildebrandlll98l . 

The uncertainties in the dust parameters used for the 
DUSTY Rts are probably the largest contribution to the 
uncertainties in the derived parameters such as Tin and 
-Mdust. With the models and the data we have, it is difficult 
to determine grain size distributions. The same SED can 
be equally well fit by very different grain populations. UV 
observations and more detailed modelling of dust emission 
are both needed to better constrain grain size distributions 
and dust chemical composition. 

6. Scaling laws revisited 

As mentioned in the Introduction, there is previous ev- 
idence that neither IR nor radio luminosity accurately 
traces the SFR in the low-luminosity low-metallicity 
galaxies studied here. In this section, we examine how well 
the canonical correlations among MIR, far-IR, and radio 
luminosities are obeyed by the BCDs in our sample. 

6.1. Radio/FIR correlations 

Several previous authors have noted that the ubiquitous 

radio to IR correlation (e.g., ICondonI 1 19921) does not 

hold for some low-meta llicitv dwarf galaxies IIDale et al.l 
I — - — 1 1 11 — - — I ■' ' ' 

|2001b; Beck ct al. 2002) or galaxies with young starbursts 

(Rousscl et al, 2003) . Such objects are characterized by 

high IR-radio ratios, and the common explanation is a 
lack of the synchrotron emission which typically domi- 
nates the radio cont i nuum of luminous la rge spirals (e.g., 
IPrice Duridri99i iRonssel et ahlEioi or free-free ab- 
sorption at 1.49 GHz ijCondon et al.lll99ll) . Both of these 
effects are at work in the BCDs in our sample, judging 
from our fitted radio indice s (see Table and previous 
detailed radi o observations (Kobulnickv & Johnson 1999; 
iTurner et aL.200Q: .Beck et al. 2002: mnt^et^aL , 20Q4a) . 

First, we examine the "g ratio" of the galaxies for 
which we have adequate data, q, defined by the logarithm 
of the ratio of the far-IR flux and the radio flux at 1.4 GHz 
llHelou et alJll98,'THCondonlll992l) . for most galaxy popu- 
lation s is very robust: (q) — 2.34 ± 0.19 (Co ndon et alJ 
Il99lh . He 2-10 and NGC 5253 have the hi ghest q val- 
ues i n our sample with q « 2.6 (see also IVacca et all 
l2002t iRoussel et all 1200.1) . These values are typical of 
ultraluminous infrared galaxies, thought to be powered 
by dense compact starbursts w hich are optically thick at 
1.49GHz ("Cond on et al]ll99ll) . The q values for Mrk33, 
IIZw40, and IIZw70 (2.5, 2.4, 2.3, respectively) are ap- 
parently representative of the normal galaxy population. 
SBS 0335—052 has a rather low q, 2.0, which is due to the 
reduced IR emission longward of 30 fim. The most "de- 
viant" galaxies (He 2— 10 and NGC 5253 on the high side, 
and SBS 0335-052 on the low side) differ > la from the 
mean value; thus our results neither confirm nor rule out 



the applicability to BCDs of the total FIR- 1.4 GHz cor- 
relations. Apparently BCDs are not extreme cases globally 
ijR.o^issel et al.lEoO^ . although radio and IR emission do 
show strong de viations from the correlation on smal l spa- 
tial scales Ce.g.. ITurner et al.ll200nl: feeck et al.ll2nn2l) . 

Next, we examine the ratios of the 60 fim and 
radio luminositi es that a re ne cessary for the appli- 
cability of the ICondonI lll992^ scaling relations for 
SFR. In order for the 60 /L(m-1.4 GHz luminosities to 
give the same SFR, th e "standard" SE D must have 
/60Mm//i.4GHz = HI (Ho pkins et ai ] l2002l) . With the ex- 
ception of SBS 0335— 052 and IIZw70, all galaxies with 
available data exceed this value by factors of 2 to 3. The 
SFRs inferred from the 60 /im luminosity are thus 2 to 
3 times higher than those from the radio, and several 
times higher (1-3 M(n vr~^) than the m edian value of 0.3 
Mq yr~^ found bv lHopkins et al.l l|2002|) for a similar sam- 
ple of BCDs. If we assume that all the radio emission is 
thermal, rather than the 10% /90% thermal/non-thermal 
mix at 1.4 GHz l)Condonll992l) . we would infer even higher 
SFRs: 2-14 Mq yr'^. TableElgives the SFRs inferred from 
the IR and radio luminosities under standard assumptions 
(see also Sect. 16.2]) . The thermal radio estimate for SFR 
is only valid for those BCDs such as NGC 5253, IIZw40, 
and SBS 0335—052 with rather flat radio spectra. The de- 
viations from the "standard" radio-IR correlation appear 
to be greater for the "active" BCDs, although our sample 
is small, and we have few passive galaxies with which to 
compare. 

These results are in contrast with those of 
iHonkins et all l|2002h who found that SFRs inferred from 
60/Ltm and 1.4 GHz luminosities are in close agreement for 
local BCDs. The main reason for the difference may be due 
to the paucity of data in their sample: of 50 sources, only 
18 were detected with IRAS at 60 ^m, and 23 at 1.4 GHz. 
Both the radio and IR data are flux-limited, with limits 
of 1-2 mJy for the radio survevs (NVS S and FIRST) and 
200 mJy for IRAS ijHopkins et al.l2002fl . The ratio of these 
limits is very close to the prediction of the IR-radio corre- 
lation, perhaps making it possible to mimic the common 
trend. Moreover, few of the lowest-metallicity galaxies are 
detected at both IR and radio wavelengths. Of the ~20 
BCDs in their sample with Z <Zq /lO, only three have 
both detections; two of these three have IR-inferred SFRs 
greater than radio ones. 

The extremely metal poor "active" BCD 
SBS 0335— 052 apparently does follow the IR-radio 
correlation; the SFRs inferred from the 60 /im and 
1.4 GHz luminosities are the sa me, 0.2 Mp,yr~^ . Ho wever, 
the true SFR is -IMgyr-i l)Hunt et alJ l2004aj) . The 
60 ^m luminosity is anomalously low because of the 
absence of cool dust, and the radio emission is low 
because of the strong thermal contribution and free-free 
absorption at 1.4 GHz. Th is is yet another aspect of the 
"conspiracy" described bv iBelll (2003) : in low- luminosity 
dwarfs, "the radio-IR correlation is linear not because 
both radio and IR emissions track SFR, but rather 
because both radio and IR emissions fail to track SFR in 
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inde pendent, but coincid entally quite similar, ways" (see 



mae penaent. Put coincia e 
also IPrice fc Dur'i3ll992|) . 



Starburst age can be another factor in determining 
the validity of the IR-radio correlations, and the abil- 
ity of these luminosities to accu rately trace SFR. In their 
studies of star-forming gal axies. iR.oiissel et al. ( 2003) and 
ICannon fc SkillmanI l|2004l) have found that the radio con- 
tinuum emission in you ng (<10Myr) starbursts ten ds to 
be pure thermal. Indeed, Cannon fc SkillmanI (|20^ have 
suggested that the observed radio spectral index be used 
as an age indicator. Heavily embedded star clusters will 
have positive spectral indices at ages ^ 1 Myr, associ- 
ated with optica lly thick thermal bremsstrahlung (e.g., 
iBeck et allbOOOl) . Then as they become increasingly opti- 
cally thin, the spectral index will flatten to —0.1. When 
supernovae (SN) and supernova remnants (SNRs) begin to 
contribute to the radio emission after ~10Myr, the radio 
index will steepen. Evolved starbursts in spiral galaxies 
in which cosmic-ray diffusion processes have taken over 
( /t 50-100 Myr) will have the "canonical" radio spectral 
index of - -0.8 ijCondo Jl flflj . 



6.2. MIR/FIR correlations 

The availability of ISO data for large numbers of galaxies 
has prompte d the use of MIR lum i nosities to trace star for- 
matio n (e.g.. lR,oussel et allEioil: iForster Schreiber et all 
l2004|) . Moreover, the relative constancy of the 7 to 15/Ltm 
ratio in normal star-forming galaxies has led to the formu- 
lation of "standard" SE Ds, and the potential application 
to hig h-redshift sources ijPale et alJl2001aHDale fc Heloul 
l2002|) . The good correlation between PAH strength and 
cold dust in normal star-forming disks ((Haas et alJl2002l) 
is one of the tenets of such formulations, but it is not 
clear that PAH strength can be used to unambiguously 
trace star formation. In Sect. we found that the "ac- 
tive" BCDs were not well fit by the starburst template, 
M82. Here we examine the implications in the context of 
MIR-diagnostics of SFRs and MIR flux ratios. 

Because of the absence of PAH features or AFEs in 
the active BCDs, and the presence of warm dust, such 
galaxies do not follow the usual trends for MIR emission. 
In NGC 5253, II Zw 40, and SBS 0335-052, the 15 to 7 /im 
flux ratio is ^ 3 times larger than would be expect ed from 
the correlations presented in lRoussel et al.l (|200ltl . These 
authors emphasize that their MIR correlations are only 
applicable in star-forming regions which are dominated 
by AFEs, but without spectral information it is difficult 
to assess this a priori. Moreover, the ratio of 15 to 850 /im 
flux in NGC 5 253. II Zw 40, and He 2- 10 exceeds that pre- 
dicted by the IPale fc HeloiJ ((2fl02f ) models by a factor of 
~ 3; in SBS 0335-052, where we have no SCUBA detec- 
tion, even the lower limit is greater than the "standard" 
prediction. We have no data for the 7 /im flux of He 2—10, 
but the remaining active BCDs show a similar trend for 
the 7 to 850 iim flux ratio. 



Again because of the suppression of AFEs in the ac- 
tive BCDs, the SFR s predicted by their 7/im luminosity 
((R,oiissel et aLlEciOll) are factors of ^ 9 lower than those 
predicted by the (thermal) radio emission. In IIZw40, the 
deficit is ^50. Their 7 /um SFRs also do not compare well 
with those inferred from the 15 /Ltm luminosities, and can 
be as much as a factor of ^ 10 lower. Table [S] gives the 
SFRs inferred from the MIR, FIR, and radio luminosities, 
and enables a comparison of the different indicators which 
hold for other samples. It is clear that the active BCDs in 
our sample are highly discrepant, relative to normal star- 
forming galaxies. On the other hand, certainly II Zw 70 
and probably Mrk 33 show SFRs that scale as usual. 

Another assessment of the validity of "standard" SEDs 
to BCDs is the 24 /im q ratio (|ADDleton et alJl2004^ . The 
24 /im-to-radio correlations in the Spitzer Extragalactic 
First Look Survey give a mean ratio at zero redshift of 
(724 = log(F24/Fi.4GHz) = 0.84 ± 0.28. For the 4 ac- 
tive BCDs in our sample, this ratio is exceeded at lev- 
els of >4cr: (924) = 2.0 ± 0.16. Mrk33 and IIZw70 
have l ower q^a, (ci9a) = 1.6 , still higher than the mean 
value. lAppleton et alJ (|2004l) find a larger dispersion in 
(724 than in the equivalent 70 /im value, and suggest that 
there is much more intrinsic dispersion in the IR/radio 
ratio at 24 /im than at longer wavelengths. The SED vari- 
ations of the kind observed in our sample, due either to 
active/passive dichotomy or to metallicity, may be causes 
of such dispersion. 

7. Age and dust masses 

In Sect. Ifi.ll w e dis cussed the hypothesis of 
ICannon fc SkillmanI (|2004|) that the shape of the ra- 
dio spectrum indicates starburst age. The IR portion of 
the spectrum also depends on the age of the starburst 
because of the time variation of dust self-absorptio n and 
the UV radiation field Ce.g.. iTakeuchi et allEool . The 
effects of metallicity are more difficult to assess, both 
because of the small size of our sample and because of 
the dearth of data for the "passive" BCDs. Nevertheless, 
we can use our sample to gain insight into possible 
dust-formation mechanisms and evolution. 

The origin of dust in galaxies is directly connected to 
their star- formation history (SFH); dust is ultimately pro- 
duced by the various byproducts of stellar evolution. In 
very young ( ^ 100 Myr) or primordial starbursts. Type 
II supernovae (SNe) are the only possible dust production 
mechanism. Later, hot massive stars (e.g., Wolf-Rayet) 
and, still later, cooler less massive stars (Asymptotic Giant 
Branch, AGBs) condense dust in their stellar winds. Dust 
production in these various phases has been estimated by 
iMorgan fc Edmunds (2003); at least 100 Myr are needed 
to achieve the conditions for dust nucleation in the stellar 
evolution of low- and intermediate-mass stars. 

We wish to examine whether the dust masses in- 
ferred from our DUSTY SED models can be attributed 
to the current starburst, and whether metallicity plays 
a role. To do so, we neglect the contribution of AGBs, 
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Table 5. SFRs" from various luminosity indicators 



Name 








SFR^ 4GHz 


SFRi 4GHz 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


IZwl8 








0.04 


0.35 


IIZw40 


0.1 


1.1 


1.2 


0.6 


5.4 


SBS 0335-052 


0.2 


1.8 


0.2 


0.2 


1.8 


IIZw70=UGC9560 






0.4 


0.3 


2.4 


NGC 5253 


0.1 


0.9 


0.7 


0.2 


1.9 


Mrk33=Haro2 






3.4 


1.6 


13.5 


He 2-10 




2.6 


2.5 


0.8 


6.9 



"Mo yr- 

b A 



'^Total SFR, assuming the relation in' 


HoDkins et aljj 


"^Total SFR, assuming the relations in 


Condon 


(199J) 



Assuming the radio emission is purely thermal free-free at 10000 K l|Condonlll99l . 



which at 100 Myr are roug hly 1% of the total dust yield 
llMorgan fc Edmundsl 12003 ) . We also neglect Wolf-Rayet 
stars ([Morgan fc Edmunds 1200,'^ . because such massive 
stars are rare, and their dust yield is low, ^ 10~^ Mo/star 
J^^^hgnk^^aOH^Q^ )■ We therefore use the models of 
IHirashita et al.1 l)2002aK based on the dust yi elds from 
type II SNe given bv iTodini fc Ferraral l)200lj) . to esti- 
mate the amount of dust generated by the current burst. 
However, we first need an estimate of SFRs. Because of 
the difficulties described in Sect. it is not clear how 
we can reliably deduce SFRs for the galaxies in our sam- 
ple. As a "best guess", we adopt the thermal radio es- 
timate from Col. (6) in Table |5l for the three galaxies 
with approximately flat radio spectral indices (IIZw40, 
SBS 0335-052, NGC 5253), and the 60/m estimate (Col. 
4) for the remaining ones. 

At ages of 3-5 Myr, with the exception of 
SBS 0335-052, Type II SNe are not capable of pro- 
ducing the amount of dust inferred from our models; 
at 5 Myr, we would expect 2 x lO'^ Mq of dust from 
the IHirashita et a l. ( 2002 3) ni o dels a nd > 10^ Mq from 
those bv iMorgan fc E dmundsl l|2003l) . These masses are 
100-1000 times lower than the 10^ Mq we find from our 
models of the active BCDs. Because of the presence of 
Wolf-Rayet stars (see Sect. 0J, the present burst in all 
galaxies except II Zw 70 must be young; but with the 
exception of SBS 0335— 052, all the galaxies for which 
we have dust mass estimates have NIR colors which 
indicate an evolved stellar population. Hence, the dust 
heated by the present burst almost certainly results 
from previous episodes of star formation. For the SFRs 
of our sample, the dust yields from SNe alone at or 
beyond ages of ^ 30 Myr are more than sufficient to 
produce the ~ 10^ Mj^ of dust inf erred from our models 
(Hirashita ct al. 2002at iMorgan fc Edmunds 2003). 

On the other hand, the low dust mass we infer for 
SBS 0335—052 can easily be ascribed to the current young 
burst. Indeed, the dust yields are rather high for the 
- 10^ Mq of dust inferred from our SED. In this BCD the 
uncertainty in the SED itself makes the dust estimate par- 
ticularly problematic. However, should our dust mass be 



correct, the burst could be extremely young ( ;S 2 Myr), 
even though the radio spectral index is slightly nonther- 
mal (~ —0.3). Alternatively, either the dust production of 
a s ingle Type II SN could be lower than that assumed 
bv IHirashita et al"l (2Q023), or other mechanisms could 
be important such as outflows ijLisenfeld fc Ferrara '^lOOS*) 
or dust destruction by SNe shocks Ce.g.. iDwe*^^ Scaloi 
^980). Obviously on the basis of one galaxy it is im- 
possible to ascertain whether metallicity, youth, or com- 
pact size/density (the "active/passive" dichotomy) are 
the primary factors in shaping the unusual SED of 
SBS 0335-052. More MIR and FIR observations of low- 
metallicity BCDs are needed to resolve the question. 

8. Metallicity and dust-to-gas ratio 

One of the foundations of models of dust formation and 
evolution in galaxies has been the dust-to-gas ratio V 
('Lisenfcld fc F erraral Il 998': 'Ed mundd lioOll: Ijames et all 

M 



2002; Hirash ita et al.ll2002 b'). Here we evaluate V for our 
sample galaxies with the aim of reassessing trends of V 
with metallicity. 

The gas mass in low-metallicity BCDs consists pri- 
marily of atomic hydrogen. Molecular gas in the form 
of CO is notoriously difficult to detect in low-metallicity 
objects, and of our sample only the highest metal- 
licity objects, NGC 5253, He 2-10, Mrk33, have been 
detected in the CO (1-0) transition jBaas et al." 

im ..Meier 

5253 has also 

been detected i n a MIR line with ISO, although 
II Zw 40 was not l|RigODOulou et al. 2002) . These CO mea- 



the CO (1 - 0) transition j^; 

Kobii]nicky elin Il99,4 iBarone et alJ l20(Ml: 
2001120021: rRravo-AlfaTO et alJl2004l) . NGC 5 



[} ' 1994 
er et alJ 



ature (Kobulnickv fc Skillman"l995'; 


'van Zee et al.lll998l: 


iThuan et al. 1999a: 


Kobulnickv fc Skillman.l99fi) enable 



us to calculate total gas masses. The molecu lar gas mass 
in NG C 5253 is 0.4% that of the atomic gas jMeier et alJ 
|2002|) . but more than twice the atomic gas mass in th e 
solar-metallicity galaxy He 2—10 ijKobulnickv et alJl995|) . 
The lack of CO detections in IIZw40 and SBS 0335-052 
could imply that the gas in these BCDs is entirely atomic, 
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dust destruction efficiency which begins to be effective 
when oxygen abundance 12 + \og{0/H) reaches ^ 8. 
Their models neglect outfl ow, and are complementary to 



7.5 8 8.5 

12 + Log(0/H) 

Fig. 7. Dust-to-gas mass ratios \ogD versus oxygen abun- 
dance 12 -f \og{0/H). The four galaxies for which we 
have DUSTY models and thus dust masses are shown 
as open stars; other data fron i the literature are also 
plotted [op en circles are fr o m Ijames et all ()2002() and 
X from .Lisenfeld fc Ferraral lIlQQsj)]. Open "s quares show 
data for IIZwTO jLisenfeld fc Ferraralll998|) and Mrk33 
l) James et al.ll2002j) . both of which are in our sample. The 
vertical dotted line indicates solar metallici ty. The solid 
and s hort-dashed lines are the models by ( James et al.l 
1)2002(1 ■ and the long-dashed and dot-dashed lines the mod- 
els for various dus t destruction efficienc ies 
bv lHirashita et al.l (l2002tJ) . 



or that CO is not a good trac er of H2 at these sub-solar 
metallicities fe.g.. lTavlor et al.i.1998'1 . The dust-to-gas ra- 
tio T) is then derived by combining our model results for 
dust mass, and reporting all quantities to our assumed 
distance. 

Figure [7| shows graphically the known correlatio n be- 
tween V and oxygen abundance l|lssa et alJ Il990f) ; our 
new determinations are shown as open stars, da ta from 
the literature as open circ les l|james et al.ll200^ and x 
l|Lisenfeld fc Ferraral Il998|) ^: Seyfert galaxies and galax- 
ies for which there are no IRAS detections are not in- 
cluded. We have also indicated with open squares the two 
galaxies in our sample which follow the standard starburst 
template, and for which we did not have sufficient data 
to pe rform a DU STY fit: IIZwTO JLisenfeld fc Ferraral 
Il998|l and Mrk33 f James et al."'2002'). NGC 5253 appears 
in the Ijames et al.. (,2002 ') sample with logI? = — 2.92, in 
very good agreement with our value logi? = —2.94. 

Model predictions for the P-metallicity correlation are 
also s hown in Fig. 171 liEdmundsl l200ll iHirashita et"al] 
l2002h(l . The IHirashita et^ l|2002h^ models, shown as 
long-dashed and dot-dashed curves, vary because of the 



iLisenfeld fc Ferraral ^^1998^ who instead attribute the vari- 
ation in T> to variati o ns in t he m ass outflow rate. B oth the 
ILisenfeld fc Ferraral l|l998(l and IHirashita et"!!] l|2002bll 
models predict a no n-linear trend for T > — O/H, unlike 
lEdmundsl l|200l[l and Ijames et all l|2002() who predict a 
linear one, proposing that the interstellar dust mass is an 
approximately constant fraction of the ISM metal abun- 
dance. A solid line indicates their prediction for dust pro- 
duction by SNe, and the short-dashed one represents that 
by only evolved low- and intermediate-mass stars. It is ap- 
parent from Fig. that none of the models reproduces the 
tren ds oi T> — O/H obse rved in our a nd other BCD sam- 
ples. |james et all l(2002|) argue that ILisenfeld fc Ferraral 
(1998D underestimated the cool-dust component because 
of the unavailability of sub-mm fluxes; this would push 
the X 's upwards in Fig. Nevertheless, our estimates of 

V which take into account the cool dust in a self-consistent 
way are equally deviant from a linear model. In addition 
to outflow and dust destruction mechanisms, intermittent 
SF H and age probably pl ay important roles in governing 

V l|Hirashita et alJl2002b|) . However, without significantly 
more data at low metallicities {12 + log{0/ H) ^ 8), it will 
be difficult to decide betwe en the linear / non-li near alter- 
natives. As pointed out bv Ijames et alJ l|200S " 



^ FoUowing lHirashita et aT ('2002b'l , we have eliminated their 
factor of 0.5 in their HI mass derivation. 



dust de- 
tections in low-metallicity star forming galaxies are badly 
needed. 

9. Implications for photometric redshifts 

The determination of the spectroscopic redshift of dis- 
tant submm/mm sources has been extremely challeng- 
ing due to their very weak and uncertain optical counter- 
parts. However, the problem of the redshift determination 
has been partially solved by using their far-IR-to-radio 
SED, comp ared to local templates, to obtain a phot omet- 
ric redshift l|Yun fc Carilll2002t 'Hughes et al.'2002'l. This 
method has been tested to successfully recover the spec- 
troscopic redshifts with a reasonable accuracy (zphot — 
Zspcc < 0.3, at 1 < z < 3). Generally, even the simple 
submm-to-radio spectral index (a^^ ^^^^) provides a goo d 
approximation of the actual redshift ( Carilli fc Yunlll999(l . 
These radio-submm photometric redshifts have been cal- 
ibrated by using templates of local, massive and chemi- 
cally evolved starburst galaxies (e.g. M82, Arp220; here 
we focus on the non-AGN cases), which are characterized 
by far-IR bumps peaking at ~60-100/im and a radio tail 
dominated by synchrotron emission. Such templates are 
certainly adequate for the submm-mm s ources at z^ 2.5 
identified by SCUBA and MAMBO ( Chapman et alJ2003l: 
IPannerb aijgr^g^^^ |2Q0 ^ , \yhich are massive and metal- 
rich fe.g.. lTecza et al.ll20oi iGenzel et al.ll2003|l . However, 
the starburst galaxies which will be detected at very high 
redshift (z > 6) by the forthcoming submm facilities (and 
ALMA in particular) will probably be characterized by 
compact star forming regions and low metallicities, sim- 
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ilar to some of the local BCDs discussed in this paper. 
As consequence, the SED of primordial galaxies at z > 6 
will likely be characterized by warmer infrared bumps and 
thermal radio emission. The use of "standard" (M82-like) 
SEDs applied to primordial galaxies with SEDs similar 
to the BCDs in our sample would yield deceptively high, 
significantly overestimated photometric redshifts (up to a 
factor of 2). 

Another related factor which is likely to affect the 
SED of starburst galaxies at z > 6 is that at such early 
epochs (age of the universe <1 Gyr) AGB stars have not 
yet evolved to produce dust. At z > 6 the only viable 
mecha nism for dust product ion is through Type II super- 
novae l|Maiolino et al.ll2004|) . SN dust has different prop- 
erties than "standard" dust (smaller grains), resulting in 
warmer SEDs. As discussed in Sect.|Hl a SN component of 
dust is also likely to contribute to the war m SED of the 
BCD SBS 0335-052 (iHirashita et a l.'2002aV Photometric 
redshifts for primordial galaxies calibrated by using the 
SEDs of local BCDs will be presented in a forthcoming 
paper. 

10. Conclusions 

We have used new and archival data to derive the optical- 
to-radio spectral energy distributions for a sample of seven 
low-metallicity BCDs. Such galaxies are regarded as the 
best local laboratories for high redshift, primordial galax- 
ies, forming stars for the first time. In particular, their 
infrared and radio emission provide important constraints 
on the early stages of star formation and on the associated 
dust properties. 

The SEDs of BCDs have revealed a wealth of infor- 
mation as well as unexpected, important properties. The 
main results can be summarized as follows: 

1. The BCD SEDs deviate significantly from the stan- 
dard templates of "classical", evolved and massive 
starburst galaxies in several ways. In particular: (i) the 
radio spectrum is generally flatter, indicative of ther- 
mal emission; (ii) the shape of the infrared bump 
generally implies warmer dust temperatures (with 
SBS 0335— 052 being an extreme case in this respect); 
(Hi) the mid-IR spectra show little or no Aromatic 
Features in Emission (e.g., PAHs). 

2. The different spectral properties of BCDs relative to 
"normal" starburst galaxies also yield strong devia- 
tions from the "scaling laws". In particular, BCDs de- 
viate substantially from the radio/FIR and MIR/FIR 
relations. One of the consequences of such deviations 
is that the SFR inferred from any of these indicators 
(radio, FIR, MIR) may be wrong by large factors. 

3. The origin of the different SEDs of BCDs with re- 
spect to "classical" starburst galaxies is not clear yet. 
Very likely the very compact and dense nature of the 
star forming regions in some BCDs (the "active" ones) 
plays a major role, and is probably responsible for 
heating the dust to higher temperatures and for the 



absence of AFEs. Age may also be important, partic- 
ularly in the dominance of the thermal emission over 
synchrotron in the radio. A young age is also probably 
related to the compactness of the star forming regions 
discussed above. The different origin of dust in the 
early stages of star formation (supernovae vs. evolved 
stars), resulting in different dust properties (typically 
smaller grains), is another factor which is likely to play 
an important role in the shape of the IR SEDs. The in- 
fluence of metallicity is not clear, and additional data 
are definitely required to further investigate its effects 
on the SEDs of star- forming galaxies. 

4. We have modelled the IR SEDs by means of a dust 
radiative transfer code, which has allowed us to de- 
rive the infrared luminosities, dust radial distribution, 
masses and temperatures. An important result is that, 
with the notable exception of SBS 0335— 052, the ob- 
served dust mass cannot be accounted for by the dust 
produced during the short episode of star formation 
traced by the still ongoing burst; most of the dust must 
have been produced by supernovae or AGB stars from 
previous episodes of star formation. Another impor- 
tant result is that the dust-to-gas ratio is not a linear 
function of the metallicity, nor does the observed rela- 
tion follow the prediction of any of the more sophisti- 
cated models. The latter result suggests that the dust 
production mechanisms in different metallicity envi- 
ronments need to be reconsidered. 

5. Finally, the significant deviations of the BCD SEDs 
with respect to "classical" evolved and massive star- 
burst galaxies prompt for a revision of the photomet- 
ric redshift techniques to be applied to the primordial, 
very high redshift galaxies, expected to be detected by 
facilities such as ALMA and SKA. 
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Appendix A: SED data and references 

The fluxes used in the SEDs together with their references 
are given in Tables Al through A7 for He 2-10, IIZw40, 
IIZw70, IZwl8, Mrk33, NGC5253, and SBS 0335-052, 
respectively. Wavelengths in the tables are given as fim up 
to the "Radio" region, where they are given in cm. The 
values reported here have not been corrected for Galactic 
extinction. The cross-identification of the data references 
follow the individual data tables. 

Data references 
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Table Al. SED for He2-10 



Region 


A (,,m ) 


Flax (mjy) 


Error (mJy) 


Ref. 


OPT+NIR 


0.36 


19.66 


1.81 


Al (61") 


OPT+NIR 


0.44 


35.76 


3.29 


Al (61") 


OPT+NIR 


0.44 


44.61 


3.70 


A2 


OPT+NIR 


0.55 


50.25 


4.63 


Al (61") 


OPT+NIR 


0.61 


112.86 


9.35 


A2 


OPT+NIR 


0.64 


63.76 


5.87 


Al (61") 


OPT+NIR 


0.79 


77.72 


7.16 


Al (61") 


OPT+NIR 


1.2 


103.78 


9.56 


Al (30") 


OPT+NIR 


1.2 


178.53 


2.30 


A2 


OPT+NIR 


1.2 


210.92 


19.42 


A3 


OPT+NIR 


1.6 


126.07 


11.61 


Al (30") 


OPT+NIR 


1.6 


183.57 


2.87 


A2 


OPT+NIR 


1.6 


227.30 


20.93 


A3 


OPT+NIR 


2.2 


104.13 


9.59 


Al (30") 


OPT+NIR 


2.2 


164.42 


3.63 


A2 


OPT+NIR 


2.2 


182.63 


16.82 


A3 


MIR 


10.8 


600.0 


40.0 


A4 


MIR 


10.8 


720.0 


95.0 


A5 


MIR 


11.7 


850.0 


20.0 


A6 


MIR 


11.7 


880.0 


9.0 


A7 


MIR 


12.0 


1093.0 


109.0 


A8 


MIR 


25.0 


6511.0 


651.0 


A8 


FIR- mm 


60 


24080.0 


2408.0 


A8 


FIR- mm 


100 


26400.0 


2640.0 


A8 


FIR- mm 


450 


580.0 


150.0 


This paper 


FIR- mm 


850 


140.0 


15.0 


This paper 


FIR- mm 


0.1 


5.9 


3.30 


A9 


FIR- mm 


0.3 


6.2 


0.60 


A9 


Radio 


2.0 


10.9 


1.00 


A7 


Radio 


2.0 


21.1 


1.20 


A9 


Radio 


2.8 


40.0 


15.00 


A9 


Radio 


3.4 


35.0 


8.00 


A9 


Radio 


36 


24.3 


2.40 


A9 


Radio 


3.6 


29.0 


3.00 


A9 


Radio 


4.8 


45.0 


10.00 


A9 


Radio 


60 


55.0 


10.00 


A9 


Radio 


62 


12.6 


4.30 


A9 


Radio 


11.1 


75.0 


35.00 


A9 


Radio 


21.4 


69.2 


6.90 


A9 


Radio 


73.5 


240.0 


80.00 


A9 


Radio 


92.2 


100.0 


20.00 


A9 



Table A3. SED for IIZwTO 



Region 


A (p,m ) Fl»: 


K (mJy) 


Error (mJy) 


Ref. 


OPT+NIR 


0.44 


4.94 


0.55 


A 19 


OPT+NIR 


0.44 


4.98 


0.60 


A2 


OPT+NIR 


0.44 


5.22 


0.48 


A20 


OPT+NIR 


0.55 


5.41 


0.65 


A2 


OPT+NIR 


0.55 


6.04 


0.56 


A20 


OPT+NIR 


0.64 


5.92 


0.71 


A20 


OPT+NIR 


0.64 


6.26 


0.58 


A 19 


OPT+NIR 


1.2 


5.81 


0.32 


A2 


OPT+NIR 


1.6 


6.03 


0.44 


A2 


OPT+NIR 


2.2 


4.98 


0.41 


A2 


FIR- mm 


25.0 


133.0 


26.6 


A2 


FIR- mm 


60 


714.0 


50.0 


A2 


FIR- mm 


100 


1237.0 


123.7 


A2 


FIR- mm 


850 


20.0 


UL 


This paper 


Radio 


2.8 


2.7 


0.20 


A 17 


Radio 


2.8 


4.1 


2.00 


A 17 


Radio 


2.8 


4.4 


1.00 


A 17 


Radio 


6.1 


3.0 


1.00 


A 17 


Radio 


6.3 


2.5 


1.10 


A 17 


Radio 


63 


3.4 


1.00 


A 17 


Radio 


20.5 


4.3 


1.00 


A 17 


Radio 


21.2 


5.0 


0.60 


A 17 


Radio 


19.1 


6.5 


0.80 


A17 


Radio 


91.7 


11.2 


2.00 


A 17 


Radio 


92.2 


12.5 


3.00 


A 17 



[Al] Johansson, I. 1987, A&A, 182, 179 

[A2] NASA Extragalactic Database (NED) 

[A3] Noeske, K. G., Papadcros, P., Cairos, L. M., & 

Pricke, K. J. 2003, A&A, 410, 481 
[A4] Telesco, C. M., Dressel, L. L., & Wolstencroft, R. D. 

1993, ApJ, 414, 120 
[A5] Vacca, W. D., Johnson, K. E., & Conti, P. S. 2002, 

AJ, 123, 772 

[A6] Sauvage, M., Thuan, T. X., & Lagage, P. O. 1997, 
A&A, 325, 98 



Table A2. SED for IIZw40 



Pegion 


A (,™ ) 


Flux (m.Jy) 


Error (mJv) 


Ret. 


OPT+NIR 


0.36 


1.91 


0.07 


AlO (166") 


OPT+NIR 


0.36 


0.47 


0.04 


All 


OPT+NIR 


0.43 


3.29 


0.09 


AlO (166") 


OPT+NIR 


0.44 


1.58 


0.15 


All 


OPT+NIR 


0.44 


2.74 


0.33 


A2 


OPT+NIR 


0.54 


6.94 


0.13 


AlO (166") 


OPT+NIR 


0.55 


4.98 


0.60 


A2 


OPT+NIR 


0.55 


5.31 


0.10 


A12 


OPT+NIR 


0.55 


5.77 


0.53 


All 


OPT+NIR 


0.64 


6.99 


0.64 


All 


OPT+NIR 


0.64 


8.56 


0.16 


A12 


OPT+NIR 


0.79 


10.06 


0.93 


All 


OPT+NIR 


0.79 


8.76 


0.16 


A12 


OPT+NIR 


0.82 


12.22 


0.34 


AlO (166") 


OPT+NIR 


1.2 


8.40 


0.31 


This paper (20") 


OPT+NIR 


1.2 


8.63 


0.44 


A2 


OPT+NIR 


1.6 


9.13 


0.34 


This paper (20") 


OPT+NIR 


1.6 


10.01 


0.68 


A2 


OPT+NIR 


2.2 


7.54 


0.65 


This paper (20") 


OPT+NIR 


2.2 


8.99 


0.33 


A2 


MIR 


7.0 


50.0 


50.0 


This paper (ISO archive) 


MIR 


8.0 


80.0 


50.0 


This paper (ISO archive) 


MIR 


10.0 


120.0 


50.0 


This paper (ISO archive) 


MIR 


10.1 


200.0 


20.0 


A13 


MIR 


10.4 


220.0 


40.0 


A13 


MIR 


11.7 


240.0 


40.0 


A13 


MIR 


12.0 


300.0 


50.0 


This paper (ISO archive) 


MIR 


12.0 


460.0 


50.0 


A14 


MIR 


13.0 


350.0 


50.0 


This paper (ISO archive) 


MIR 


14.0 


400.0 


50.0 


This paper (ISO archive) 


MIR 


15.0 


530.0 


50.0 


This paper (ISO archive) 


MIR 


20.0 


1000.0 


200.0 


A13 


MIR 


25.0 


1910.0 


200.0 


A14 


FIR-mm 


60 


6610.0 


700.0 


A14 


FIR-mm 


100 


5800.0 


900.0 


A14 


FIR-mm 


450 


1550.0 


350.0 


This paper 


FIR-mm 


850 


50.0 


10.0 


This paper (40") 


FIR-mm 


850 


90.0 


10.0 


This paper (total) 


FIR-mm 


1300 


13.3 


3.1 


A15 


Radio 


1.2 


18.0 


4.00 


A16 


Radio 


2.0 


12.0 


3.00 


A16 


Radio 


2.8 


20.0 


1.00 


A16 


Radio 


6.0 


22.0 


0.10 


A16 


Radio 


20.1 


30.5 


1.50 


A17 


Radio 


20.7 


30.0 


0.50 


A16 


Radio 


21.4 


34.2 


1.40 


A18 


Radio 


92.2 


38.0 


4.00 


A17 
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Region 


A („m ) 


Flux (m.Jy) 


Error (mJy) 


Ret. 


OPT+NIR 


0.439 


1.10 


0.05 


A21 (20") 


OPT+NIR 


0.44 


1.62 


0.06 


AlO 


OPT+NIR 


0.555 


1.19 


0.05 


A21 (20") 


OPT+NIR 


0.64 


0.98 


0.06 


AlO 


OPT+NIR 


0814 


0.76 


0.03 


A21 (20") 


OPT+NIR 


1.2 


0.76 


0.03 


A21 (20") 


OPT+NIR 


1.6 


0.64 


0.02 


A21 (20") 


OPT+NIR 


2.2 


0.51 


0.02 


A21 (20") 


FIR-mm 


850 


2.5 


UL 


This paper 


FIR-mm 


1200 


5.5 


UL 


F. Bertoldi, priv. comm. 


Radio 


2.8 


2.1 


0.60 


A16 


Radio 


3.5 


0.6 


0.04 


A22 


Radio 


3.5 


0.8 


0.09 


A22 


Radio 


6.0 


1.6 


0.20 


A23 


Radio 


6.2 


1.0 


0.18 


A22 


Radio 


6.2 


0.5 


0.07 


A22 


Radio 


6.3 


1.6 


0.20 


A16 


Radio 


21.0 


1.2 


0.19 


A22 


Radio 


21.0 


1.6 


0.15 


A22 


Radio 


21.0 


2.6 


0.90 


A16 



[A7] Beck, S. C, Turner, J. L., & Gorjian, V. 2001, AJ, 

122, 1365 
[A8] ISSA/OASIS 

[A9] Kobulnicky, H. A. & Johnson, K. E. 1999, ApJ, 527, 
154 

[AlO] Marlowe, A. T., Meurer, G. R., Heckman, T. M., & 

Schommer, R. 1997, ApJS, 112 
[All] Heisler, C. A. & Vader, J. P. 1994, AJ, 107, 35 
[A12] Telles, E. & Terlevich, R. 1997, MNRAS, 286, 183 



20 Hunt, Bianchi, & Maiolino: SEDs of low-metallicity BCDs 



Table A5. SED for Mrk33 








Region 


A (,im ) Flu 


^ (mjy) 


Error (mJy) 


Ref. 


OPT+NIR 


0.44 


16.96 


1.56 


A2 


OPT+NIR 


0.44 


18.60 


3.42 


A24 


OPT+NIR 


0.44 


18.77 


0.69 


A 19 


OPT+NIR 


0.55 


13.46 


0.62 


A25 


OPT+NIR 


0.55 


17.42 


O.SO 


A 25 


OPT+NIR 


0.64 


21.91 


2.22 


A 19 


OPT+NIR 


0.64 


28.09 


2.59 


A24 


OPT+NIR 


1.2 


43.66 


1.21 


A2 


OPT+NIR 


1.2 


47.44 


4.37 


This pBp6r 


OPT+NIR 


1.6 


52.55 


4.84 


This paper 


OPT+NIR 


1.6 


53.53 


1.73 


A2 


OPT+NIR 


2.2 


39.59 


3.65 


This paper 


OPT+NIR 


2.2 


44.22 


1.63 


A2 


MIR 


12.0 


210.0 


31.0 


A2 


MIR 


25.0 


1050.0 


30.0 


A2 


FIR-mm 


60 


4770.0 


42.0 


A2 


FIR-mm 


100 


5990.0 


129.0 


A2 


FIR-mm 


450 


170.0 


UL 


This paper 


FIR-mm 


850 


50.0 


10.0 


This paper 


Radio 


1.2 


10.0 


3.00 


A 16 


Radio 


1.3 


3.2 


1.80 


A 16 


Radio 


2.0 


16.0 


0.06 


A26 


Radio 


2.8 


7.1 


0.90 


A 16 


Radio 


36 


4.9 


0.05 


A26 


Radio 


6.0 


11.0 


0.05 


A26 


Radio 


6.1 


31.0 


6.00 


A 16 


Radio 


63 


12.0 


1.00 


A 16 


Radio 


11.1 


24.0 


6.00 


A 16 


Radio 


20.0 


20.0 


0.08 


A26 


Table A7. SED for SBS 0335-052 






Region 


A (pm ) Fin 


X (mJy) 


Error (mJy) 


Ref. 


OPT+NIR 


0.36 


0.58 


0.02 


A35 


OPT+NIR 


0.44 


0.70 


0.01 


A35 


OPT+NIR 


0.55 


0.80 


0.01 


A36 


OPT+NIR 


0.64 


0.73 


0.01 


A35 


OPT+NIR 


0.79 


0.45 


0.00 


A36 


OPT+NIR 


1.2 


050 


0.02 


A37 


OPT+NIR 


1.6 


0.37 


0.02 


A37 


OPT+NIR 


2.2 


0.44 


0.02 


A37 


OPT+NIR 


38 


058 


0.11 


A38 


MIR 


5.3 


0.9 


0.1 


A39 


MIR 


60 


2.9 


0.1 


A39 


MIR 


67 


5.5 


0.4 


A40 


MIR 


70 


5.8 


0.1 


A39 


MIR 


77 


8.9 


0.5 


A40 


MIR 


81 


9.8 


0.1 


A39 


MIR 




11.6 


0.1 


A39 


MIR 




11.3 


0.1 


A39 


MIR 


10.2 


12.0 


0.1 


A39 


MIR 


108 


10.3 


5.1 


A41 


MIR 


109 


16.0 


0.1 


A39 


MIR 


11.3 


20.7 


0.8 


A40 


MIR 


11.3 


18.4 


0.1 


A39 


MIR 


11.8 


24.0 


0.1 


A39 


MIR 


12.0 


20.5 


0.6 


A40 


MIR 


12.0 


25.1 


0.1 


A39 


MIR 


13.3 


37.3 


0.1 


A39 


MIR 


14.3 


41.7 


0.1 


A39 


MIR 


14.9 


35.2 


0.9 


A40 


MIR 


15.1 


42.8 


0.1 


A39 


MIR 


170 


45.3 


0.1 


A39 


MIR 


18.0 


49.3 


0.1 


A39 


MIR 


19.0 


55.3 


0.1 


A39 


MIR 


20.1 


55.7 


0.1 


A39 


MIR 


21.0 


50.2 


15.1 


A41 


MIR 


21.0 


57.5 


0.1 


A39 


MIR 


21.6 


60.1 


0.1 


A39 


MIR 


22.0 


70.0 


11.0 


A39 


MIR 


22.5 


59.0 


0.1 


A39 


MIR 


23.1 


59.0 


0.1 


A39 


MIR 


24.6 


60.8 


0.1 


A39 


MIR 


25.1 


61.9 


0.1 


A39 


MIR 


25.9 


63.0 


0.1 


A39 


MIR 


26.6 


64.2 


0.1 


A39 


MIR 


27.4 


64.2 


0.1 


A39 


MIR 


28.4 


65.0 


0.1 


A39 


MIR 


29.4 


64.6 


0.1 


A39 


FIR-mm 


30.2 


62.6 


0.1 


A39 


FIR-mm 


31.6 


62.0 


0.1 


A39 


FIR-mm 


65 


112.0 


21.0 


A41 


FIR-mm 


65 


44.0 


20.0 


This paper 


FIR-mm 


850 


5.0 


UL 


This paper 


Radio 


0.130914 


11.00 


UL 


A42 


Radio 


0.267672 


2.30 


UL 


A42 


Radio 


1.3 


0.5 


0.08 


A43 


Radio 


2.0 


0.4 


0.07 


A43 


Radio 


3.5 


0.7 


0.04 


A43 


Radio 


6.2 


0.8 


0.09 


A43 


Radio 


20.1 


0.4 


0.13 


A42 


Radio 


20.5 


0.5 


0.06 


A43 



Table A6. SED for NGC 5253 



Region 


A (,™ ) 


Flnx (mJy) 


Error (mJy) 


Ret 


OPT+NIR 


0.36 


74.76 


6.88 


All 




OPT+NIR 


0.44 


129.84 


11.96 


All 




OPT+NIR 


0.44 


191.17 


21.12 


A2 




OPT+NIR 


0.55 


174.22 


16.04 


All 




OPT+NIR 


0.55 


242.72 


26.82 


A2 




OPT+NIR 


0.64 


273.24 


25.16 


All 




OPT+NIR 


1.2 


147.27 


13.56 


A27 


(48") 


OPT+NIR 


1.2 


381.70 


8.44 


A28 




OPT+NIR 


1.6 


163.15 


15.02 


A27 


(48") 


OPT+NIR 


1.6 


412.86 


10.64 


A28 




OPT+NIR 


2.2 


133.53 


18.44 


A27 


(48") 


OPT+NIR 


2.2 


329.58 


12.75 


A28 




OPT+NIR 


3.8 


145.01 


13.35 


A27 


(r.'5) 


MIR 


4.5 


125.0 


30.0 


A29 




MIR 


7.8 


865.7 


233.7 


A30 


(8'.'2) 


MIR 


8.3 


1076.20 


99.10 


A27 


(r.'5) 


MIR 


8.6 


1153.9 


115.4 


A30 


(«'2) 

(r.'5) 


MIR 


9.4 


1215.50 


111.90 


A27 


MIR 


9.6 


1065.9 


149.2 


A30 


{g!2) 


MIR 


10.0 


1297.7 


77.9 


A30 


(8'.'2) 


MIR 


10.3 


1200.0 


20.0 


A27 


(7.'5) 


MIR 


10.3 


1372.00 


126.30 


A27 


(7.'5) 


MIR 


10.4 


1539.0 


169.3 


A30 


(S'2) 


MIR 


10.8 


1340.0 


40.0 


A4 




MIR 


11.4 


1573.9 


188.9 


A30 


(8'.'2) 


MIR 


11.7 


2200.0 


140.0 


A31 




MIR 


12.0 


1960.0 


300.0 


A29 




MIR 


12.0 


2720.0 


33.0 


A27 


(7.'5) 


MIR 


12 


2215.60 


204.00 


A2 




MIR 


12.4 


2041.8 


449.2 


A30 


(8'.'2) 


MIR 


14.9 


3900.0 


400.0 


A29 




MIR 


17.4 


7910.20 


728.40 


A27 


(7'.'5) 


MIR 


18.7 


2900.0 


190.0 


A31 




MIR 


20.0 


6098.7 


914.8 


A30 


(8'.'2) 


MIR 


25.0 


13060.0 


48.0 


A2 




FIR-mm 


60 


30000.0 


54.0 


A2 




FIR-mm 


100 


30920.0 


95.0 


A2 




FIR-mm 


450 


2200.0 


UL 


This 




FIR-mm 


850 


180.0 


20.0 


This 




FIR-mm 


1200 


114.0 


4.0 


A29 




Radio 


2.0 


54.0 


5.00 


A32 




Radio 


3.6 


58.0 


10.00 


A33 




Radio 


6.0 


75.0 


10.00 


A2 




Radio 


6.0 


49.0 


5.00 


A32 




Radio 


6.2 


90.0 


12.00 


A2 




Radio 


20.0 


55.0 


10.00 


A32 




Radio 


21.4 


90.0 


40.00 


A32 




Radio 


73.5 


128.0 


10.00 


A34 
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Appendix B: Dust mass from DUSTY 

Here we set out the equations used to infer the dust mass 
from the output given by DUSTY. Assuming that the dust 
is distributed in a spherical shell with a power-law radial 
density profile p(r) = Cir~P, the dust mass can be written 



as: 



M, 



dust 



AttCi 



3-p P ^ 3 
In^ p = 3 



(B.l) 



where Ci is a constant, and Rin and i?out are the inner 
and outer shell radii. 

DUSTY provides the value of r^, the optical depth at 
frequency ly, which can be expressed as: 

T„ = I X„{r)dr 



-Rout /•a.2 

dr / da na^ 

Rin Jai 



(B.2) 



where i refers to the grain chemical composition type. X^, 
is the absorption coefficient per unit length; Qi^v{a) the 
extinction efficiency factor (ratio of optical to geometrical 
cross section) for type i grains with size (radius) a; ni{a, r) 
the number density distribution in radius for grains of type 
i at radius r of size (radius) a; and 7ra^ the geometrical 
cross section of a grain with size (radius) a. ai and 02 
are the minimum and maximum grain sizes, respectively. 
ni[a,r) can be written as two functions which separate 
the mass distribution radial dependence from the grain 
size: 71^(0, r) = Ni(r)T{a). Assuming the grain sizes are 
distributed as a power law, !F{a) = Aa~'^, where A is 
such that J^^ J-'{a) da — 1. Ni(r) is the number density 
of grain type i at radius r. 

We can now write pi the mass density at radius r of 
grain type i as: 



/4 

Si ( gTra ) ni{a,r) da 



= -TTS,N,{r)AB 



(B.3) 



where B = J^^ a^ ' da and Si is the mass density of the 
material composing grains of type i. But Pi{r) can also be 
expressed as Pi{r) = fip{r) = fiCir~P, where fi is the 
fractional mass density abundance for grains of chemical 
type i. We can then solve for Ni{r) in terms of the constant 
Ci: 

- CiD,r-P (B.4) 

where Di = 4^^^^^ , and depends only on the grain prop- 
erties. Eg nation IB . 21 can then be rewritten in terms of Di 
and Ci: 

fRoMt ra,2 

dr^N,{r) I daTTa^Q^^^{a)T{a) (B.5) 



drCiY^Di 
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da (7ra2)Q^^^(a)Aa-« (B.6) Bressan, A., Silva, L., & Granato, G. L. 2002, A&A, 392, 
ai 377 

Cannon, J. M. & Skillman, E. D. 2004, ApJ, 610, 772 



/a.2 
-1 



T^CiA{ i-p f^-K Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 

\\n^ p=l] 345 245 

CariUi, C.L. & Yun, M.S. 1999, ApJ, 513, L13 

(B.7) Calzetti, D., Meurer, G. R., Bohlin, R. C, Garnett, D. R., 
Kinney, A. L., Leitherer, C, & Storchi-Bergmann, T. 

Here {Qi i,{a)) is a mean value, since it has been taken out ^^^"^^ AJ, 114, 1834 

of the integral Chapman, S. C, Windhorst, R., Odewahn, S., Yan, H., & 

^, .'. n n ■ ^ Consehce, C. 2003, ApJ, 599, 92 

The quantities i?i„, i?out, P, ai, 02, g are mput; hence ' ' ^ ' ' 

we can calculate the constants A and B. DUSTY requires ^^^'^^ ^ ^1^^^' ^001, ApJ, 556, 562 

the specification of the fractional number abundances of S-, Lacey, C. G., Baugh, C. M., & Frenk, C. S. 2000, 

the grains Fi = NJY^.Ni, rather than the fractional MNRAS, 319, 168 

mass density abundances However, the two quantities ^^^don, J. J. 1992, ARA&A, 30, 575 

are related through: J' Anderson, M. L., & Helou, G. 1991, ApJ, 

^ 376, 95 

5,N, _ S,F, Conti, R S. 1991, ApJ, 377, 115 

h - Pi/P- S^N, ~ 5^F^ ^ ' Cox, A. L., Sparke, L. S., Watson, A. M., & van Moorsel, 

G. 2001, AJ, 121, 692 

Therefore, with the density 5, for each grain type i, Dale, D. A., Helou, G., Contursi, A., Silbermann, N. A., 

A can be calculated. Assigning a mean extinction effi- & Kolhatkar, S. 2001a, ApJ, 549, 215 

ciency factor {Q,) to each grain type, we can finally re- Dale, D. A., Helou, G., Neugebauer, G., Soifer, B. T., 

late to the constant Ci, necessary to derive the dust Prayer, D. T., & Condon, J. J. 2001b, AJ, 122, 1736 

mass as given in Equation |BT] In our calculations, we Dale, D. A. & Helou, G. 2002, ApJ, 576, 159 

have assumed the following physical propert ies for the Dannerbauer, H., Lehnert, M. D., Lutz, D., Tacconi, L., 

grains: 5^^,r' = 1.81 gcm-^ llKim et alJll994l): (5siiicate = Bertoldi, F., Carilh, C, Genzel, R., & Menten, K. M. 

3.3gcm-3 llDraine fc Le3ll984l): 6„..„w,,. = 2.26gcm-3 2004, ApJ, 606, 664 

llDraine fc Leelll984l): and (Q,) = 2 for aU grain types de Jong, R. S. 1996, A&A, 313, 377 

(see|Spitzei|[l978|). de Grijs, R., et al. 2004, MNRAS, 352, 263 

Desert, F.-X., Boulanger, F., & Puget, J. L. 1990, A&A, 
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